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D.  L  Graf,  C.  R.  Biyth,  and  R.  S.  Stemmler 


ABSTRACT 

The  diffraction  effects  expectable  from  CaMg  car- 
bonates having  a  variety  of  abnormal  sequences  of  basal 
cation  population  have  been  calculated  by  programming  the 
University  of  Illinois  digital  computer  ILLIAC  I  for  the  most 
general  of  the  Hendricks -Teller  equations.  Cases  involv- 
ing two  and  three  kinds  of  multiple  layers  (packets)  with 
correlation  extending  to  the  immediately  preceding  packet 
(s  =  range  of  correlation  =  1)  aretreatedas  Markov- 1  proces- 
ses. The  random  packet  succession  (s=0)  is  a  particular 
solution  of  the  (s=l)  case.  In  a  second  paper,  to  be  pub- 
lished later,  the  two-packet  (s=2)  case  is  formulated  in 
terms  of  a  two-packet  Markov-2  model,  which  is  then  trans- 
formed  into  a  four-packet  Markov- 1  model  but  not   solved. 

The  two-packet  (s=l)  model  reproduces  the  basal 
diffraction  features  observed  for  synthetic  protodolomite 
powders  formed  near  room  temperature.  Three-packet  (s=  1) 
models,  used  as  a  device  for  introducing  more  complex  cor- 
relations, are  required  for  categorizing  the  replacive  dolo- 
mite rhombohedra  found  in  sedimentary  limestones. 

The  unambiguous  evidence  for  CaMg  carbonate  mixed 
layering  available  thus  far  is  contained  in  a  few  zero-level 
a-axis  Weissenberg  diagrams .  In  the  best  defined  example, 
the  nonideality  of  the  basal  cation  plane  succession  appar- 
ently involves  the  distribution  of  pure  Mg  layers  and  Ca- 
substituted  Mg  layers,  rather  than  the  insertion  of  extra 
Ca  layers. 


INTRODUCTION 

Variations  in  the  composition  and  crystal  structure  of  the  carbonate  min- 
erals comprising  limestone  and  dolomite  are  known  to  affect  the  physical  and  chem- 
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ical  properties  of  these  rocks  in  industrial  utilization.    As  part  of  a  continuing  long- 
range  program  of  limestone  and  dolomite  studies,  we  examine  here  the  X-ray  dif- 
fraction effects  expectable  for  carbonate  crystals  having  one  particular  kind  of 
crystallographic  defect. 

One-dimensional  disorder  has  been  reported  for  a  number  of  carbonate  com- 
pounds.   Coplanar  carbonate  ions  in  alternate  (010)  carbonate  layers  of  UO2CO0 
(orthorhombic)  point  in  opposite  directions  in  the  ideal  structure,  but  diffuse  streak- 
ing in  the  direction  of  the  b*  axis  in  single-crystal  X-ray  diffraction  diagrams  of 
some  crystals  indicates  occasional  regions  in  which  anions  point  the  same  way  in 
successive  layers  (Christ,   Clark,   and  Evans,    1955).    This  is  an  example  of  Jagod- 
zinski'  s  (1949)  "Lagenfehlordnung,  "  in  which  lattice  planes  of  the  same  composi- 
tion are  disoriented  in  their  position  with  respect  to  each  other.     Parisite,   2CeFC03' 
CaC03  (rhombohedral),   normally  is  made  up  of  basal  CeFC03    and  CaC03  layers 
in  regular  Ce:Ce:Ca  sequence,  but  diffuse  reflections  drawn  out  parallel  to  c*  in- 
dicated disordered  layer  sequences  in  two  crystals  examined  (Donnay  and  Donnay, 
1953).     Planar  carbonate  groups  must,   from  optical  data,   lie  nearly  perpendicular 
to  the  layers  in  this  compound. 

Basal  reflections  with  asymmetrical  and  even  bimodal  intensity  distribu- 
tions, as  contrasted  with  essentially  symmetrical  profiles  for  reflections  with  no 
e-axis  component,  have  been  reported  in  powder  and  single-crystal  X-ray  diagrams 
of  rhombohedral  carbonates  near  dolomite,  CaMg(C03)2,   in  composition  (Graf  and 
Goldsmith,    1956;  Goldsmith  and  Graf,    1958;  Goldsmith  and  Heard,    1961).     In  the 
ideal  dolomite  structure,  coplanar  carbonate  ions  lie  in  basal  anion  planes  inter- 
leaved with  basal  cation  planes  alternately  populated  by  the  two  cations.     Inten- 
sity relations  for  the  aberrant  dolomites  appear  on  casual  observation  to  be  normal 
not  only  among  those  reflections  with  no  basal  component  but  also  among  those  of 
various  orientations  whose  intensity  results  principally  from  oxygen  scattering. 
Thus,   there  is  no  evidence  for  significant  departure  of  anion    positions  or  orienta- 
tions from  those  characteristic  of  dolomite,  and  localized  in-basal-plane  cation 
ordering  of  the  huntite  type  (Graf  and  Bradley,    1962),  which  requires  a  different 
anion  arrangement  from  those  of  calcite  or  dolomite,  becomes  improbable.    The  do- 
lomite basal -re  flection  anomalies  were  attributed  to  a  defective  succession  of  ca- 
tion-plane populations  (Jagodzinski'  s  (1949)  "Artenfehlordnung") .    Goldsmith  and  Graf 
(1958)  pointed  out  that  the  c:a  ratios  for  calcite,  dolomite,  and  magnesite  differ 
enough  so  that  an  appropriate  mixture  of  solid  solutions  will  yield  basal  reflec- 
tions attenuated  with  respect  to  those  with  no  basal  component.    However,  they 
cited  single -crystal  diagrams  in  which  the  difference  in  reflection  profiles  was  too 
great  to  be  explained  in  this  way. 

Before  further  diffraction  study  of  these  dolomites,  criteria  must  be  estab- 
lished for  defining  the  nature  of  cation-plane  successions  (layer  disorder)  and  cat- 
ion-plane  populations    (point  disorder).     In  seeking  these  criteria,  we  assume 
throughout  that  solid  solution  within  a  cation  plane  takes  place  at  random  cation  sites, 

Two  methods  are  available  for  attempting  to  define  the  nature  of  mixed-lay- 
er successions.     One -dimensional  Fourier  transforms  of  the  observed  diffraction 
effects  (MacEwan,    1953,    1956)  give  a  direct  statement  of  the  proportions  and  thick- 
nesses of  the  contributing  layer  types  and  some  idea  of  the  randomness  of  the  lay- 
er succession.    In  general,  the  greater  the  number  of  orders  that  can  be  measured, 
the  more  useful  the  results  given  by  such  transforms.     Some  difficulty  is  encoun- 
tered in  choosing  baselines  and  recognizing  spurious  peaks  in  the  computed  trans- 
forms, but  what  is  more  important,  no  intermediate  orders  of  significant  intensity 
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can  be  omitted,  whether  a  single  measurement  of  maximum  intensity  and  corres- 
ponding cZ-spacing  is  made  for  each  order  or  a  series  of  such  measurements  to  take 
account  of  asymmetric  profiles.     These  transforms  are  clearly  suited  to  the  anal- 
ysis of  basal  reflection  data  obtained  from  single-crystal  diffractometry. 

The  obvious  disadvantage  of  the  second  available  procedure,  the  calcula- 
tion of  theoretical  basal-reflection  intensity  distributions  for  a  series  of  mixed- 
layer  models,   is  the  large  number  of  computed  curves  often  needed  to  test  all  rea- 
sonable possibilities.    However,  because  of  our  interest  in  fine-grained  synthet- 
ic and  natural  CaMg  carbonate  powders,   for  which  only  the  |00-  6|  and    |00-  12  j 
profiles  from  the  full  sequence    {oO-s},     {oO-g},     [oO-9},     {oo-izj,    {oo-is}... 
can  be  resolved,  this  second  approach  is  obviously  necessary.    The    |00'3|  re- 
flection is  too  weak  to  be  observed  readily  in  powder  X-ray  diagrams,  even  of  ideal 
crystals;    |00"9J    is  masked  by  the  much  stronger  |ll'6  |  reflection,  and  higher 
orders  are  either  too  weak  or  masked.    The  weak  |  00*  i.  \  reflections  with  Ji  odd 
are  order  reflections  that  owe  their  existence  to  the  segregation  of  Ca  and  Mg  in- 
to alternate  basal  cation  layers.    Zero  intensity  of  these  reflections  resulting  from 
randomizing  the  cation-plane  succession  can  be  inferred  if  stronger,  nonbasal  or- 
der reflections  such  as   |01-5|  and  |02*l|  cannot  be  detected. 

It  will  be  seen  that  as  the  amount  of  mixed  layering  becomes  less,   higher 
orders  of  basal  reflection  must  be  examined  in  order  to  describe  it,   and  such  ex- 
amination beyond  |00-12|  is  possible  only  with  single  crystals. 

Hendricks  and  Teller  (1942)  developed  mixed-layer  matrix  equations  for 
several  models  of  increasing  complexity.    We  have  programmed  the  University  of 
Illinois  digital  computer  ILLIAC  I  (Graf,   Blyth,  and  Stemmler,    1957,    1958)  for  the 
most  elaborate  of  these  equations,  the  repeated  solution  of  which  would  be  quite 
impractical  by  hand  calculation.    The  two  fundamental  repeat  distances  involved 
in  considering  the  CaMg  carbonates  are  essentially  the    loO*  12j-    interplanar  spac- 
ings  for  MgC03  and  CaC03,  the  distances  between  adjoining  basal  Mg  and  car- 
bonate planes  and  between  adjoining  basal  Ca  and  carbonate  planes.    We  have 
used  1.253  and  1.422  A,  respectively,  for  these  two  repeat  distances,   but  several 
sources  of  error  in  the  choice  should  be  noted: 

(1)  The  preferred  value  for  MgCOg,   measured  after  these  calculations  were 
well  advanced,   is  1.251  A  (Graf,    1961). 

(2)  The  variation  of  the  unit  cells  of  materials  intermediate  between  calcite 
and  magnesite  is  not  perfectly  linear  (Goldsmith,  Graf,   and  Heard,    1961). 
Thus,  d  /qo.i2\    ^°^  cation-ordered  material  of  dolomite  composition  is 
1.334  A,   that  for  essentially  cation-disordered  material  is  1.337  A,  and 
the  mean  of  the  1.253  and  1.422  A  values  used  in  these  calculations, 

1 .3375  A,  approximates  a  mean  of  the  actual  1 .23  A  and  1 .45  A  values 
for  calcium  and  magnesium  packets  in  dolomite. 

(3)  Very  fine-grained  carbonate  precipitates  formed  near  room  temperature  have 
enlarged  unit  cells  tentatively  ascribed  to  hydration,  e.g.  the  sedimentary 
magnesite  described  by  Graf,   Eardley,  and  Shimp  (1961)  with  d  |oo-12|  ~ 
1.268  A.  ' 

The  result  of  the  two  repeat  distances  being  so  nearly  alike  is  that  carbon- 
ate mixed  layering  is  often  expressed  as  the  asymmetry  or  broadening  of  reflections 
that  do  not  shift  appreciably.    Evaluation  of  experimentally  obtained  profiles  thus 
inevitably  involves  other  factors  such  as  particle  size  and  residual  strain  that  af- 
fect peak  width  and  asymmetry.    We  have  not  attempted  a  rigorous  evaluation  of 
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experimental  measurements.    Instead,   for  comparison  with  the  results  of  calcula- 
tions, we  use  already  published  material  and  a  few  tentative  analyses  of  experi- 
mental data  that  assume  the  effects  of  point  and  layer  disorder  on  Kctj  peak  shapes 
to  be  dominant  over  those  from  other  causes . 

The  Hendricks -Teller  intensities  become  infinite  at  discrete  maxima  and 
are  presumably  distorted  to  a  lesser  degree  at  angles  close  to  maxima  and  for  lay- 
er successions  only  slightly  removed  from  perfect  regularity.     In  general,  there- 
fore, we  place  prime  reliance  upon  the  computed  positions  of  maxima  rather  than 
upon  the   predicted  relative  intensities  at  these  points.    We  make  but  little  use  of 
predicted  mixed-layer  broadening  as  a  diagnostic  parameter  because  of  the  diffi- 
culty of  separating  it  experimentally  from  broadening  resulting  from  particle  size, 
residual  strain,  and  the  presence  of  a  range  of  particle  compositions.    There  is  an 
obvious  need  for  a  parameter  reflecting  asymmetry,  and  we  attempt  to  use  centroids 
measured  from  Hendricks -Teller  intensity  distributions,  realizing  the  limitations 
imposed  by  distortion  of  intensities  at  sharp  maxima.    Evaluation  of  bulk  compo- 
sition is  made  by  examining  reflections  with  no  c-axis  component. 

The  behavior  of  reflections  with  no  basal  component  merits  some  discus- 
sion.    For  a  powder  made  from  a  crystal  consisting  of  magnesite  grown  onto  an 
(00 -1)  calcite  substrate,   one  would  observe  the  complete  set  of  diffraction  maxima 
of  calcite  and  that  of  magnesite.     In  dolomite,  which  can  be  considered  an  alter- 
nating sequence  of  basal  one -cation-plane  units  of  calcite  and  magnesite,   there 
has  been  contraction  of  the  lateral  repeat  distance  of  the  calcite  units,  with  smal- 
ler percentagewise  basal  expansion,   and  concomitant  vertical  contraction  and  lat- 
eral expansion  of  the  magnesite  units.    The  flexibility  of  in-plane  orientation  of 
carbonate  groups  in  dolomite  that  permits  this  rearrangement  is  a  function  of  the 
oxygen  "thermal"  ellipsoid.     One  set  of  sharp  diffraction  maxima  is  observed.     The 
extent  to  which  lateral  adjustment  operates  in  more  irregular  successions  to  leave 
the  reflections  with  no  basal  component  undistorted  is  not  known.     Obviously,  there 
must  be  some  small  upper  limit  to  the  number  of  Ca  (or  Mg)  basal  planes  that  can 
locally  exist  as  nearest  cation  plane  neighbors  in  dolomite,  without  an  intervening 
cation  plane  of  the  other  kind,  and  still  leave  such  reflection  profiles  symmetrical. 
There  is  a  second  unknown  limit  beyond  which  these  reflections  will  show  two  max- 
ima.    These  two  limits  for  reflections  with  no  basal  component  will  not,   in  general, 
be  the  same  as  for  the  basal  reflections,   because  profiles  of  the  latter  respond  to 
differences  in  scattering  power  of  successive  basal  cation  planes  as  well  as  dif- 
ferences in  spacing. 

The  group -theoretical  approach  of  Jagodzinski  (1954)  would  be  well  suited 
to  treating  three-dimensional  models  that  could  predict  profiles  of  various  nonbas- 
al  reflections,   but  at  present  we  do  not  have  sufficient  experimental  data  to  justi- 
fy such  an  analysis.    We  anticipate  that  it  would  be  most  useful  in  describing  high- 
ly segregated  interlayered  successions,  with  all  reflections  asymmetrical  or  dou- 
bled, that  would  ordinarily  be  passed  off  as  mixtures  of  two  or  more  solid  solutions. 

C.  R.  Blyth,   of  the  Department  of  Mathematics,   University  of  Illinois,   par- 
ticularized for  the  two-  and  three-packet  cases  the  general  development  contained 
in  Hendricks'   and  Teller's  Section  4.     R.  S.  Stemmler,   at  present  a  postdoctoral 
research  associate  at  the  National  Bureau  of  Standards  under  National  Academy  of 
Sciences  —  National  Research  Council  sponsorship,   programmed  the  calculations 
for  ILLIAC  I.    A  brief  discussion  of  this  programming  is  given  in  the  Appendix. 

We  are  particularly  indebted  to  the  Graduate  College  of  the  University  of 
Illinois  for  generous  use  of  ILLIAC  I  and  to  W.  F.  Bradley  for  many  helpful  and 
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stimulating  discussions.     Irmgarde  Martin,   Sharon  Wolfe,  John  Dyni,  and  Roland 
Stemmler  have  at  various  times  shared  with  us  the  tedium  of  processing  large  mas- 
ses of  numbers, 

Donald  L.  Graf  was  employed  by  the  Illinois  State  Geological  Survey  when 
most  of  the  work  reported  in  this  paper  was  carried  out,   but  has  since  moved  to  the 
Department  of  Geology  and  Geophysics,   University  of  Minnesota  at  Minneapolis. 
A  first  draft  of  the  manuscript  was  written  by  him  during  the  tenure  of  a  research 
fellowship  extended  by  the  Committee  on  Experimental  Geology  and  Geophysics  of 
Harvard  University.     During  the  preparation  of  the  final  manuscript  at  the  Univer- 
sity of  Minnesota,  he  was  supported  by  National  Science  Foundation  grant  GA-453. 


A  TWO-PACKET  MARKOV- 1  MODEL 

The  derivation  of  the  two-packet  Markov- 1  model  is  given  in  detail  in  order 
to  show  the  dependence  of  the  various  quantities  involved.    This  model  is  merely 
Hendricks'   and  Teller's  general  solution  particularized  for  two  kinds  of  multiple 
layers . 

We  use  a  formula  for  diffracted  intensity  (James,    1950,  chap.  2;  Buerger, 
1950,  chap.  3),  which  assumes  that  the  crystal,  although  containing  a  large  num- 
ber of  unit  cells,   is  small  enough  so  that  there  is  no  appreciable  rescattering  or 
absorption  of  radiation.     For  a  layer  of  atoms  of  one  kind,  e.g.  Ca,   impinged  upon 
by  X-radiation  that  is   polarized  with  the   electric  vector  perpendicular  to  the  plane 
of  incidence, 

V^^^fQane^X/mc^sin^  (1) 

where  n  is  the  number  of  Ca  atoms  per  unit  area,  X  is  the  wavelength  of  the  inci- 
dent radiation,  9  is  the  angle  of  incidence,   f^^  is  the  atomic  scattering  factor  for 
Ca,   e  and  m  are  the  charge  and  mass  of  the  electron,  c  is  the  velocity  of  light  in 
free  space,   and  Vq^  is  the  ratio  of  reflected  to  incident  amplitude,  the  latter  as- 
sumed to  be  of  unit  value  and  zero  phase. 

Similarly,   for  a  layer  of  coplanar  carbonate  groups, 

VQQ^={fQ+  3fQ)ne^X/mc^sin  9  .  (2) 

The  number  of  cations  and  carbonate  groups  per  unit  area  in  the  respective 
basal  planes  is  the  same,   so  that  n  can  be  set  equal  to  unity  in  both  (1)  and  (2). 
Indeed,   because  relative  rather  than  absolute  magnitudes  are  desired,  the  products 
of  constants  in  (1)  and  (2)  may  be  set  equal  to  unity.    Although  X  is  a  constant  in 
any  given  calculation  or  experiment,   it  is  retained  in  the  equations  so  that  results 
will  be  expressed  in  terms  of  (sin  9)/\,  rather  than  sin^,   and  thus  will  be  appli- 
cable to  experimental  measurements  made  with  X-radiation  of  any  wavelength.    The 
phase  lag  between  amplitudes  contributed  by  planes  a  distance  a  apart  is  47ra(sin 
9)/\,  and  thus  the  resultant  amplitude  Vq^qoo  fro'"  ^  layer  of  Ca  atoms  plus  the 
carbonate  layer  immediately  preceding,   summed  at  the  Ca  layer,   is 

VcaC03=    ^[fca+    (fc  -  3fo)  exp  (-i477-a;isin  ^  )/x]   /(sin  ^  ),  (3) 
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where  a^  is  the  separation  of  adjacent  Ca  and  carbonate  planes.     Similarly, 


^MgC03~  ^ 


^Mg  "^   ^^C  "*"  '^^O^  ®^P  (-i47ra2Sin^  )/X 


'^{sm9).  (4) 


where  a2  is  the  separation  of  adjacent  Mg  and  carbonate  planes. 

Equations  (3)  and  (4),  of  course,   can  also  be  written  in  polar  form: 

VqqqOo  ~  ^Ca^^^^"^  ^    +  fQQ„X/sin0     [cos  47ra2(sin  ^  )/X -i  sin  47ra-]^(sin  ^  )/X  J 

and 

VMgC03  =  fivig^/sin  ^   +  f^Qo^^sin  0    [cos  4"'a2(sin  5  )/X-i  sin  47ra2(sin  ^  )/Xj  . 

In  an  alternate  packet  definition,  the  packets  may  be  considered  as  made  up  of 
three  layers,   CO3/2  — Ca  — CO3/2  (each  a^  apart)  and  CO3/2— Mg  — CO3/2  (each 
32  apart).     In  this  case,   summing  at  the  cation  plane, 

^CaC03^  ^Ca^^^^"  ^    '''  ^COs  ^  [  cos  4  77-a  j(sin  ^  )/Xl /sin  9 

and  similarly  for  Vn^gQQ   .     For  this  packet  definition,  which  leads  ultimately  to 
calculated  intensities  identical  with  those  obtained  for  the  two-layer  packets,  the 
amplitudes  V  are  no  longer  complex.     However,  the  imaginary  components  of  scat- 
tering do  not  cancel  for  asymmetrical  packets  used  in  the  subsequent  three-packet 
model  and,   in  any  event,  the  quantity  0  appears  in  the  calculations  in  complex 
form.     For  this  reason,  the  more  general  packet  constructions  involving  complex 
notation  were   used  throughout  the  paper. 

Because  cation  and  anion  planes  alternate  in  the  rhombohedral  carbonates, 
it  has  been  possible  to  reduce  the  problem  from  one  involving  three  kinds  of  layers 
(Ca,  CO3,  and  Mg)  to  a  two-layer  case.    We  shall  henceforth  refer  to  multiple  lay- 
ers like  CaC03  as  packets. 

Consider  a  hypothetical  series  of  CaMg  carbonates,  ranging  from  calcite 
to  magnesite  in  composition,  made  up  of  CaC03  and  MgC03  packets  without  solid 
solution  of  one  cation  in  packets  appropriate  to  the  other.    Let  f^  '  and  f^  '  be  the 
frequencies  of  the  two  kinds  of  packets.    We  adopt  the  Hendricks -Teller  model  so 
that  the  packets  form  a  stationary  Markov  chain  with  two  possible  states,   initial 
probabilities  f'   ',    f'^),  and  matrix  of  transition  probabilities 

P21P22/      ^21^        ^22^ 

where  Pg^.  is  the  probability  that  a  type  s  packet  will  be  followed  by  a  type  t  packet. 
This  model  is  Markov- 1  in  nature  — the  probabilities  are  dependent  upon  the  kind  of 
packet  immediately  preceding,  but  not  upon  still  earlier  ones.    The  reason  for  writ- 
ing the  second  expression  for  P  above  is  that  if  the  packets  were  independently  as- 
sorted,  P  would  have  this  value  with  all  Cg^=  1. 
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Because  the  matrix  P  has  row  sums  1, 


and  stationarity  requires 


Cj^f        +0^2^      =1.  (6) 


f^\/'K   f^\/'K,(^\ 


which  for  f      >0  reduces  to 

(2) 

When  f       >0,   subtraction  of  (7)  from  (6)  gives  C|2  =  C2i.     In  the  nontrivial  case 

f(l)  >0,   f(2)  >o,  p  can  be  written  in  the  form 

'"  *      cf"'   l-cf»' 

Here  c  is  a  parameter  that  can  range  over  the  open  interval    |o,   min(l/f^     , 

l/f(2))j  .    The  Dnobability  of  a  randomly  selected  pair  of  adjacent  packets  being 

different  is  2cf'     f^^'  compared  with  2f'l)f^^'  if  the  packets  were  independent.     Thus 

c  measures  the  degree  to  which  the  two  kinds  of  packets  tend  to  alternate  in  the 

crystal.     The  value  c=  1  indicates  random  assortment.     Smaller  values  of  c  indicate 

a  tendency  of  like  packets  to  cluster  together;  larger  values  of  c  indicate  a  tendency 

of  the  two  kinds  of  packets  to  alternate . 

Following  Hendricks  and  Teller  (1942),  we  write  Qi  for  the  matrix    Ipst  exp 
(-i<^(st))|  ^  ^^  L 

/ll-cf(2^exp(-i<^(ll))         cf(2)  exp  (-i^(12^ 
""^^y     cf^^^xp  (-10(21))     (i-cf(l))exp(-i0(22))y. 

where  (f}^^^>  is  the  phase  shift  between  adjacent  packets  of  types  s  and  t,   i.e. 

(/)  (^^^  =  87ra^  sin  9  j>  ^^^^^STra^  sin   9 


^  (12)^^(21)^4^(^^^   g^)  g.^     Q 


Rewriting  Qj^  as 

^21   ''22/' 


Q^=[''^^''^M  (8) 
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the  characteristic  equation  belonging  to  (8)  is 


^11-^   ^^12 


'21 


^22-^ 


=  0 


(9) 


and  the  solutions 


^=^^'Kl^   ^22V^^11^^22)'^   ^^12^21} 


(10) 


of  (9)  are  Q^'-'  and  O'    ',  the  eigenvalues  of  Q^.     In  obtaining  the  square  root  of 
the  complex  number  in  (10),  we  use  the  formula  ya  +  bi  =  x  +  iy,  where 


J  /2       T  1/2 

x=  ±  K/a    +    b    +    a 


and 


y=± 


a    +    b 


U/2 


and  X  and  y  have  like  sign  if  b>0,   unlike  if  b<0.    We  are  free  to  choose  the  order 
of  the  characteristic  roots,  and  thus  let  Q'^'  be  the  root  having  a  real  part>0,    q'     , 
the  one  having  a  real  part<0. 

Introducing  the  matrix  O, 


0= 


^11   ^12 


^21   ^22 


we  write  the  matrix  equation  OQ   =Q0 


^1   ^2 


^21   ^22 


^11   ^^12 


^21   ^522 


from  which 


^11^11  ^^2^21  =  QJJ|  ^11 

^1^12  +  ^2^22=^^        ^2 

(2) 
^2l''ll  +  ^22^21  =  ^(2)^21 

^21^12  +  ^22^22  =  ^        ^22. 


(11) 
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Each  bracketed  set  of  two  equations  in  (11)  is  a  homogeneous  system  in  two  un- 
knowns with  a  zero  determinant  because  its  Q^^'  is  a  root  of  (9).     One  a^j  in  each 
pair  of  equations  therefore  can  be  chosen  arbitrarily  (but  not  zero),   and  we  take 


a   „  =  1  and  a      =  1 .    We  then  have 


"11=^22-^ 


1) 


^22=^11-^ 


(2) 


-q 


12 


-q 


21 


(12) 


The  matrix  R  is  defined  as 


^u'n'    /v'"V"  v'"V^'\  /f"'o 


R=  OVFO"  ^  =  I 


'22 


-a 


12 


^1^22-^   ^1^22"^ 


21      22/     \  /    \  /     Vll^22-^    ^1^22-^/' 

(D*  (1) 

where  V  is  the  complex  conjugate  .of  V      .    The  diagonal  elements  of  R  are 

°h'22-' 


r(22)=.j(")U(1)V"+  a,„v'2'V'U  a,,£<2'U")*vC'+   a^v'^'V^'} 


22 


11 


22 


^1^22-^ 

Substituting  into  the  Hendricks -Teller  expression  (1942,  equation  30,   p. 
159)  for  the  average  intensity  of  scattered  radiation  per  layer,  we  have 


i^^W)  =  t 


(1) 


V 


(1) 


+  f 


(2) 


,(2) 


2       R^^V^)       ^(22)^(2)        r(11)*q(1)*       r(22)*q(2)^ 
^      I-q(^)  1-q(2)  i-o(l)*     ^        i-q(2)*     ' 


which  may  be  rearranged  to  give 

2 


W^)=f 


(I) 


V 


(1) 


,,.,     ,..2        2W^^^0^^V2 

1-Q 


Q 


(1) 


'm^''h 


1) 


+   2K(R^^^^Q^^V2 


0 


(2) 


W^^b 


(13) 


1-Q 


(2) 
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Figure  1  -  Basal  reflection  sequences  for  crystals  consisting  of  equal  numbers  of 
CaCOj  and  MgCOj  packets,  at  various  values  of  the  parameter  c,  computed 
using  the  two-packet  model. 
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where  Jpi denotes  "real  part  of."    For  c=0  and  for  c  =  Cj^ajj  ,  I^y  becomes  infinite 
at  the  (sin  Q  )/\  value  of  each  order  of  basal  reflection.     Equation  (13)  is  for  mono- 
chromatic radiation,   and  it  does  not  include  the  customary  Lorentz  or  polarization 
corrections.     To  account  for  the  electric  vector  lying  in  the  plane  of  incidence,  I^y 
should  be  multiplied  by  cos     26;  for  unpolarized  incident  radiation,   the  common 
experimental  situation,   it  should  be  multiplied  by  (1  +  cos^  29)/2.    The  pertinent 
Lorentz  geometry  is  that  of  the  single  crystal  case.    We  do  not  make  these  correc- 
tions for  the  experimentally  measured  intensities  cited  later  in  this  paper,  because 
these  measurements  extend  only  over  limited  angular  ranges. 


TWO-PACKET  MARKOV- 1  CALCULATIONS 

Figure  1  shows  the  first  four  orders  of  basal  reflection  computed  for  a  crys- 
tal made  up  of  equal  numbers  of  CaCO^  and  MgC03  packets,  for  several  values  of 
c.    The  atomic  scattering  factors  used  for  C,   O,   Mg,  and  Ca  in  this  and  subsequent 
calculations  are  those  of  Berghius  et  al.   (1955),  all  self-consistent  field  data  with 
exchange.    In  accord  with  the  empirical  observation  that  X-ray  scattering  by  these 
carbonates  is  more  nearly  that  of  neutral  atoms  than  of  ions,  the  scattering  curve 
of  Berghius  et  al.  for  Mg"*""*"    has  been  modified  at  low  values  of  (sin^)/X  toward 
that  expectable  for  neutral  Mg. 

Similar  curves  for  mixtures  of  CaC03  and  MgCOj  packets  in  the  ratio  of 
75  to  25  are  shown  in  figure  2.    This  is  a  bulk  composition  typical  of  the  more  high- 
ly magnesian  skeletal  calcites.    The  Mg  in  these  metastable  materials  has  usually 
been  considered  as  randomly  distributed  among  the  total  cation  positions,   but  if  it 
were  to  segregate  completely  into  discrete  MgC03  packets,  the  curves  of  figure  2 
show  the  diffraction  effects  expectable  for  various  distributions  of  these  packets 
among  the  CaC03  packets.    We  did  not  calculate  order  reflection  intensities  for 
this   bulk  composition  with  the  two-packet  Hendricks -Teller  program.    However, 
the  Ca2CaMg  block,   part  of  the  illustration  in  the  subsequent  section  on  three- 
packet  models,   gives  hand-calculated  intensities  for  this  composition  and  the 
maximum  possible  c,    1.33.    For  this  hypothetical  crystal,   basal  reflections  would 
occur  at  a  (sin  9)/\  interval  half  that  of  figure  1 . 

We  next  assume  that  substitution  of  Ca  into  MgC03  packets,  and  vice  ver- 
sa, can  be  treated  by  allowing  both  the  thicknesses  of  the  packets  and  the  scat- 
tering powers  of  the  cation  planes  to  vary  linearly  between  the  values  for  the  un- 
substituted  cases  and  by  changing  f^  '  and  f^^^  as  needed  to  secure  desired  bulk 
compositions.     Figure  3,  a  plot  derived  from  these  calculations,  gives  the  separa- 
tion of  the  two  maxima  in   |00-6|-  and  |  00 -121   reflections  having  a  bimodal  inten- 
sity distribution.     For  a  given  definition  of  packet  types,  variation  in  bulk  compo- 
sition has  only  a  moderate  effect  on   A  (sin  ^)/X  ;  before  a  curve  of  very  small  or 
very  large  c  has  diverged  greatly  from  the  others,   it  ends.    Although  the  graph  is 
incomplete,  it  can  be  seen  that  the  range  of  c  values  over  which  doubled  reflec- 
tions exist  becomes  smaller  as  substitution  within  individual  cation  planes  increas- 
es.   It  also  should  be  remembered  that  the  range  of  bulk  compositions  that  can  be 
treated  decreases  for  the  models  with  cation  substitution,   because  the  composi- 
tions of  the  two  packet  types  are  the  bulk  composition  limits. 

The  value  of  plots  like  that  of  figure  3  is  limited  by  the  fact  that  the  same 
A  (sin  9  )/X  results  from  many  different  possible  packet  definitions.     Figure  4  gives 
graphs  of  a  different  type  summarizing  the  positions  of  intensity  maxima  and  cen- 
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troids  of  I  00-  6  j  and   (oo-  12|  .     In  these  diagrams,  the  diagonal  straight  line  is 
the  locus  of  intensity  maxima  corresponding  to  the  maximum  possible  c  for  each 
bulk  composition.    The  maximum  possible  c  is  continuously  variable  along  the  di- 
agonal line,  reaching  a  maximum  of  2  at  the  middle  and  declining  to  zero  at  each 
end.     The  range  of  (sin  Q  )/\  values  corresponding  to  packet  successions  varying 
from  random  to  ideally  alternating  (cj^^j^  to  c=  1)  is  much  greater  for  the  higher  or- 
der   loo-  12 1  reflection,  a  point  of  considerable  practical  importance  in  trying  to 
characterize  slightly  disordered  dolomites. 

In  the  diagrams  of  figure  4,  the  centroid  position  for  a  given  bulk  composi- 
tion approaches  that  of  the  intensity  maximum  as  c  increases,   and  the  diagonal  line, 
therefore,   is  also  the  locus  of  centroid  positions  for  maximum  possible  values  of 


75%  CoCOj  packets 
25%  MgC03  packets 
Maximum    C  =    1.33 
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Figure  2  -  Basal  reflection  sequences  for  crystals  consisting  of  7  5  percent  CaC03 
packets  and  25  percent  MgCOj  packets,   at  various  values  of  the  parameter  c, 
computed  using  the  two-packet  model. 
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c.     Reflections  of  infinite  intensity  at  c^^^  thus  need  not  be  treated  explicitly  and 
offer  no  difficulty  in  constructing  the  diagrams . 

However,  the  continuous  background,   not  apparent  at  the  scale  of  figures 
1  and  2,   increases  with  packet  randomization  and  presents  a  serious  problem  in 
calculating  centroids.    At  c=  1.95  and  f^'-' =  f^'^' ,  the  point  intensity  predicted  at 
the    [OO-  IZJ    maximum  is  (in  arbitrary  units)  660,000  compared  with  30  at  the  point 
of  lowest  background  between    |oO-6|  and    |00-12J  .    At  c=  1,  the  comparable  val- 
ues are  17,  000  and  600.     Prediction  by  the  model  of  this  rising  background  into 
which  the  attenuated  odd -order  basal  reflections  finally  disappear  is  in  accord  with 
experimental  observations  in  other  systems.    Although  one  could  include  the  long 
reflection  tails  in  model  calculations  of  centroid  positions,  experimental  measure- 
ments typically  involve  pseudocentroids  in  limited  (sin  9  )/X  regions  close  to  re- 
flections,  using  an  arbitrary  nonzero  baseline.    A  similar  procedure  was  followed 
in  constructing  the  diagrams  of  figure  4.     For  each  reflection  a  section  of  the  basal 
diffraction  profile  extending  for  .060  to  .075  in  (sin  9)/\  was  used,   if  available. 

Centroid  positions  are  quite  sensitive  to  baseline  height,   i.e.  to  the  length 
of  profile  considered,  and  can  be  shifted  erroneously  if  intensity  values  are  select- 
ed at  too  great  an  interval  relative  to  the  widths  of  sharp  maxima.    As  we  said  ear- 
lier,  individual  intensities  at  sharp  maxima  are  exaggerated  and  therefore  will  domi- 
nate the  profile  with  regard  to  centroid  position.     Much  of  the  difficulty  we  have 
had  in  contouring  centroid  positions  undoubtedly  stems  from  these  factors.    Some 
of  the  centroid  lines  in  the  diagrams  of  figure  4  are  sketched,   and  even  the  best 
of  them  are  much  more  poorly  located  than  those  for  maxima. 

The  curve  of  centroid  positions  for  c=  0  in  each  diagram  is  particularly  im- 
portant because,  together  with  the  diagonal  of  the  figure,   it  delimits  the  complete 
family  of  centroid  curves.    However,  the  c=  0  case  is  that  of  epitactic  overgrowth, 
and  it  is  particularly  clear  in  this  limiting  situation  involving  two  distinct  phases 
that  intensities  must  be  corrected  for  mass  absorption  if  the  calculated  centroids 
are  to  be  comparable  with  those  of  real  crystals.     Solids  in  the  interval  from  c^^^ 
to  c=  1  are  clearly  single  phase,  requiring  no  mass  absorption  correction,  and  those 
with  increasing  segregation  as  c  becomes  progressively  less  than  1  presumably  re- 
quire larger  and  larger  corrections.    The  centroid  curves  in  figure  4,   based   upon 
uncorrected  Hendricks -Teller  intensities,  therefore  must  be  progressively  more  in 
error  for  values  of  c  less  than  1.    It  would  be  necessary  to  transform  experiment- 
ally measured  intensities  if  the  diagrams  of  figure  4  were  to  be  used  in  their  pres- 
ent form. 

The  centroid  curves  of  figure  4,  within  the  limits  of  accuracy  of  our  proce- 
dure,  intersect  at  a  common  point,   as  do  those  curves  of  maximum  intensity  posi- 
tion that  run  the  full  length  of  the  diagram.     The  centroid -curve  intersections  lie 
close  to,   but  not  precisely  upon,  the  diagonals  of  the  figures.     Likewise,   the  inter- 
section of  the  maximum-intensity-position  curves  in  each  diagram  lies  close  to  the 
intersection  of  the  diagonal  with  the  vertical  line  representing  the  midpoint  of  the 
compositional  range  under  consideration.    As  the  two-packet  types  become  more 
nearly  alike   in  composition,   the  centroid  common  intersection  point  moves  toward 
the  edge  of  the  diagram  and  the  envelope  of  centroid  curves  lying  against  the  di- 
agonal becomes  shallower.    Likewise,  the  common  intersection  point  of  the  maxi- 
mum-intensity-position curves  moves  toward  the  intersection  of  the  center  vertical 
with  the  diagonal,  and  the  range  of  c  values  included  within  the  envelope   of  con- 
tinuous curves  of  this  type  extends  to  lower  and  lower  values.    We  have  not  drawn 
common  intersection  points  for    |00'  6|   centroid  curves  on  diagrams  involving  im- 
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Figure  3  -  Relation  between  the  parameter  c  and  the  separation  of 
maxima  of  bimodal    {00-6V    and    {00-12V    intensity  distributions. 
Packet  types  are  defined  in  terms  of  the  atomic  fractions  of  Ca 
and  Mg  in  basal  cation  planes.    The  numbers  within  the  dia- 
grams give  the  Ca  atomic  fraction  of  the  total  cation  population 
of  the  bulk  crystal.    The  precise  points  of  disappearance  of  sub- 
sidiary maxima  have  not  been  determined;  individual  curves  are 
drawn  to  the  last  value  of  c  at  which  two  maxima  were  observed 
in  the  ILLIAC  calculations  available. 
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pure  cation  layers,   because  we  do  not  have  enough  computed  profiles  at  very  CaC03- 
rich  compositions  to  locate  these  points. 

Packet  thicknesses  are  clearly  not  independent  of  the  proportions   of  the  two 
kinds  of  packets,   and  we  now  seek  to  estimate  the  effect  of  this  dependence  upon 
positions  of  intensity  maxima.     Steinfink  and  Sans  (1959)  obtained  zq=  .2440  ± 
.00017  and  zq=  .2435  ±  .00031  (hexagonal  cell)  for  a  dolomite  with  a ^=  16. 119  A. 
An  earlier  determination  by  Bradley,  Burst,  and  Graf  (1953)  on  a  dolomite  with  Cq  = 
16.  0  A  was  z=  0.  243 .    A  e^  value  of  16.010  A  was  proposed  for  dolomite  of  ideal 
composition  and  structure  by  Goldsmith  and  Graf  (1958).     The  1953  z  parameter  leads 
to  Mg  — 003=  1.223  A  and  Ca— C03=  I.444  A,   and  the  1959  zq  and  zq  values  lead 
to  Mg— 003=  1.235  A  and  Ca  — 003=  1.433  A.    We  use  the  latter  set  of  values,  with 
the  reservation  that  an  uncertainty  in  z  parameter  as  a  function  of  composition  exists. 

For  packet  successions  intermediate  between  perfect  alternation  and  perfect 
segregation,  the  thickness  of  a  specific  packet  in  the  succession  will  be  a  function 
of  the  kinds  of  packets  occupying  the  first  few  positions  on  either  side.    We  as- 
sume that  the   effect  of  the  immediately  neighboring  packet  is  dominant,   so  that  av- 
erage packet  thickness  can  be  defined  as  varying  linearly  with  the  Markov- 1  para- 
meter c.     Figure  5  gives  the  results  of  such  calculations  for  a  model  involving  sin- 
gle-cation packets  and  bulk  compositions  near  dolomite.    At  c=2,  Ca — 003=  1.433 
A  and  Mg— C03=  1.235  A;   at  c=  0,  Ca  — 003=  1.422  A  and  Mg— 003=  1.253  A,   etc. 
This  model  shows  that  treating  packet  thickness  as  a  variable  produces  significant 
but  not  drastic  changes  in  the  positions  of  maxima,  especially  at  intermediate  val- 
ues of  c.     The  common  intersections  of  curves  of  maximum  intensity  position  noted 
in  figure  4  have  been  destroyed. 

Synthetic  Protodolomite  as  an  Example 

Graf  and  Goldsmith  (1956)  described  fine-grained  dolomitelike  precipitates 
formed  below  100°  C  that  had  no  ordering  reflections,  greatly  attenuated    JOO-  12  j 
maxima,  and   |00-6|   peaks  that  were  displaced  and  attenuated.    If  the  possibility 
of  bound  water  or  hydroxide  is  ignored,  these  precipitates  appear  to  correspond  to 
the  two-packet  curves  of  figure  1  with  c  between  0.5  and  1.    Thus  they  either  are 
random  packet  successions  or  an  assemblage  of  calcite  and  magnesite  slabs,  rath- 
er than  a  dolomite  with  imperfect  order.     However,  given  c/a  ratios  in  the  CaMg 
carbonates  that  range  from  3.24  to  3.42,  there  is  a  question  whether  or  not  a  mix- 
ture  of  solid  solutions  could  produce  the  diffraction  features  observed.    We  have, 
therefore,   chosen  one  sample  for  further  examination,   applying  to  it  an  empirical 
graphical  analysis. 

Figure  4  A-H  -  A  series  of  diagrams,   derived  from  two-packet  model  calculations, 
which  relate  the  positions  of    loO-ej    and    |oO-12|    centroids  and  intensity  max- 
ima to  various  values  of  c,   packet  composition,   and  bulk  composition.    Solid 
contour  lines  define  positions  of  intensity  maxima,   and  dashed  lines  give  cen- 
troid  positions.    Segments  of  solid  contour  lines  with  an  arrow  at  each  end  run 
the  full  length  of  the  diagram,   but  are  shown  only  in  part.     Each  solid  line  that 
ends  within  the  diagram  has  been  drawn  only  to  the  last  point  at  which  the  com- 
puted basal  reflection  profile  still  shows  that  particular  intensity  maximum  in 
the  ILLIAC  computations  available.     Precise  bracketing  of  the  points  of  disap- 
pearance of  the  auxiliary  maxima  would  require  considerable  additional  com- 
putation. 
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A  protodolomite  formed  by  holding  a  CaMg  carbonate  gel  of  approximately 
equimolar  composition  for  276  hours  at  100°  C  (Graf  and  Goldsmith,    1956,  run  53, 
table  3)  consists  of  spherulites  approximately  3  microns  in  diameter.    Splinters  of 
the  spherulites  are  length-fast,    so  that  the  fibers  must  be  elongated  along  the  c?- 
axis  or  at  least  along  a  direction  closer  to  c  than  to  a.     This  material  shows  no  or- 
der reflections.     The  profile  across    |00'6|  ,   obtained  by  point-counting  at  0.025° 
intervals  in  2^,   followed  by  graphical  resolution  of  the  Kct|  component  using  the 
method  of  Rachinger  (1948)  and  Pease  (1948),   is  given  in  figure  6  together  with  the 
profile  of  the  nearby    Ill-OJ-  reflection,  which  has  no  c-axis  component. 

The  protodolomite    |ll'  OJ-   reflection  in  figure  6  has  a  width  at  half  maxi- 
mum intensity  of  about  0.45°  in  2^,  compared  with  0.20°  for    {H'OJ   of  the  well 
crystallized  dolomite  from  Binnental,   Switzerland,  run  under  the  same  experiment- 
al conditions.    We  estimate,   accordingly,  that  the  effective  particle  size  across 
the  fibers  is  about  2  50  A.     This  corresponds  to  about  375  fibers  making  up  the  cir- 
cumference of  an  equatorial  section  of  a  3-micron  spherulite.    A  minimum  number 
for  comparison  can  be  obtained  from  electron  micrographs  of  the  spherulites,  which 
show  jagged  edges;  there  are  about  100  separate  projections  around  a  typical  spher- 
ulite . 

The    |00-6[  reflection,  which  is  generated  in  the  direction  of  elongation  of 
the  fibers  and  thus  presumably  involves  a  larger  crystallite  size  than    |ll'  OJ,   is 
nevertheless  seen  in    figure  6  to  have  two  maxima  and  to  be  more  attenuated  than 

The   |ir  0|    and  |00-  6|  Ka-]^  maxima  each  can  be  divided  into  a  pair  of  rough- 
ly symmetrical  profiles  of  about  the  same  angular  width,  one  of  which  has  its  maxi- 
mum essentially  at  the  IQ  value  for  CaMg(C03)2.    We  make  this  graphical  resolu- 
tion,  assuming  for  the  moment  that  the  sample  is  made  up  of  two  phases  having  these 
compositions.     The  difference  in  intensity  of  the  two  reflections  can  be  ignored  be- 
cause in  each  case  we  are  interested  in  subprofiles  within  the  reflection. 

The  CL^  versus  composition  plot  of  Goldsmith,  Graf,  and  Heard  (1961)  for 
cation-disordered  solid  solutions  is  next  used  to  estimate  the  compositions  cor- 
responding to  the  two    jll-  o|    subprofile  maxima.    By  referring  these  two  composi- 
tions to  the  Cq  plot  of  Goldsmith,  Graf,   and  Heard  (1961),  we  shift  the  two   111-  o| 
subprofiles  to  corresponding    |00-6l   subprofile  positions  and  thus  construct  a  syn- 
thetic   |00-6|    Kci;^  profile.      The  instrumental  correction  mentioned  in  the  legend  of 
figure  6  is  considered  during  these  manipulations. 

It  c^n  be  seen  from  figure  6  that  the  secondary  intensity  contribution  in  the 
synthetic   |00-  6|    profile  is  at  the  same  angular  position  as  that  in  the  observed 
|00-6|  Kflj  profile,   but  is  quantitatively  inadequate  to  explain  the  latter.     The  sec- 
ondary intensity  contributions  are,  respectively,   about  15  and  30  percent  of  the  to- 
tal Krtj  areas.     The  two-subprofile  construction  that  we  have  made  brings  the  largest 
possible  percentage  of  Kct-j^  intensity  (15  percent)  to  the  position  of  the  high  angle 
|00*6|    maximum;  the  use  of  a  larger  number  of  subprofiles  would  lower  this  per- 
centage. /         T  r  1 

In  figure  6  we  make  a  graphical  comparison  of    |03*0j   and   |00'12|  ,  which 
lie  at  nearly  the  same  angular  position.     The  reliability  of  the  Kcc^  profile  obtained 
from  the  latter  reflection  is  lessened  by  its  extreme  attenuation,  but  the  presence 
of  a  subsidiary  higher  angle  component  is  still  clearly  indicated  by  a  break  in  slope . 
Dividing    |03-0|  Kcc,  into  a  subprofile  at  the  angular  position  for  CaMg(C03)2  and 
a  residual  subprofile  at  the  lower  angle,  to  be  consistent  with  our  treatment  of    |ll-  o| 
Kdi,   gives  subareas  of  28  percent  and  72  percent,    substantially  different  than  the 


28 


ILLINOIS    STATE    GEOLOGICAL    SURVEY    CIRCULAR    408 


I  mol  %  CoCOj  substitution 

29  in  degrees,  CuKS  rodiotion 


Graphical  resolution  of  the  Kct^   and  Ka2  components  of 


Figure  6  - 

mite.  Individual  measured  intensities  are  shown  by  filled  cir- 
files .  The  constructed  |00-6J  profile  has  for  convenience  been 
position  of  the  Ktf  Intensity  maximum  of  cation-ordered  dolomite 
in  Graf  (1961),  and  A'  is  the  position  of  the  related  Ka^  inten- 
imum  of  largely  cation-disordered  material  of  the  composition 
and  Heard  (1961).  An  instrumental  correction  for  the  diffractom- 
has  been  incorporated  in  locating  A,  A'  ,  and  A". 
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(ll-oj,    (oo-e),     {o3-o},  and  {oO-12}   for  a  synthetic  protodolo- 
cles,  and  the  dashed  lines  are  empirical  subdivisions  of  the  Knj  pro- 
plotted  at  reduced  amplitude.     In  each  profile,  A  is  the  computed 
having  the  composition  CaMg(COo)2.   based  upon  the  unit  cell 
sity  maximum.    A"  is  the  computed  position  of  the  Kffj  intensity  max- 
CaMg{C03)9,   obtained   from  the  unit  cell  given  in  Goldsmith,  Graf, 
eter  of  0.09     in  IQ ,  determined  from  measurement  of  a  Si  standard. 
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15-85  percent  split  for    |l  1-  o|    Kct^.     The  synthetic    |oO-  12J    profile  obtained  on 
continuing  the  analysis  has  too  much  diffracted  intensity  at  the  high-angle  end  of 
the  profile  and  none  in  the  middle. 

By  assuming  the  sample  to  consist  of  a  mixture  of  three  or  more  solid  solu- 
tions and,   if  necessary,   assuming  the  solid  solutions  nearer  dolomite  in  composi- 
tion to  have  lesser  half  widths,   the  observed    |00-12J    profile  could  be  matched 
more  closely  and  the  relative  importance  of  the  subprofiles  within  |  03-  o\   and    111-  o\ 
brought  into  better  agreement.     In  the  process,  however,  the  disparity  between  con- 
structed and  observed    |00-6|   profiles  would  be  increased. 

A  simple  mixture  of  solid  solutions  appears  inadequate  to  explain  the  dif- 
fraction features  exhibited  by  this  sample,  although  we  have  not  considered  pheno- 
mena connected  with  fibrous  aggregates,   such  as  differential  absorption.    The  count- 
er assumption  is  that  the  asymmetry  of    |ll-  OJ    and    |03-  o|    can  be  ignored  and  the 
diffraction  features  of    |00-6|   and   1 00 -121    attributed  solely  to  mixed  layering. 
The  (sin  ^  )/\  separations  within  the  two  latter  reflections,   if  we  assume  the  break 
in  slope  of  the   |00-12|  Ka^  profile  to  correspond  to  a  subsidiary  maximum,  are 
about  .0030  and  .0056,   respectively.     Referring  to  the  two  diagrams  of  figure  3, 
and  estimating  in  each  the  position  of  a  diagonal  line  delimiting  the  field  in  which 
curves  exist,  we  see  that  the  value  of  c  for  the  protodolomite  sample  must  be  be- 
tween 0  and  .  1  and  that  there  must  be  40  to  45  mol  percent  substitution  in  both 
kinds  of  cation  layers.    The  sample  is,  therefore,  not  much  different  than  a  mixture 
of  two  solid  solutions  near  dolomite  in  composition. 


MARKOV- 1  CASES  INVOLVING  THREE  KINDS  OF  PACKETS 

The  (00- X)  series  intensity  distributions  seen  in  single -crystal  X-ray  dia- 
grams of  some  natural  dolomite  crystals  (Goldsmith  and  Graf,    1958,   pi.  2)  differ 
from  any  of  those  predicted  by  the  two-packet  model.    The  (00-3)  and  (00-9)  reflec- 
tions have  nonzero  intensity  for  crystals  with  asymmetrical  or  bimodal  (00-  6)  and 
(00-  12)  reflections,  and  one  (at  least)  of  the  odd -order  reflections  (00-  9)  is  clearly 
asymmetrical  or  bimodal  in  the  various  diagrams. 

The  presence  of  part  of  the  diffracted  (00-6)  and  (00  -12)  intensity  at  (sin 
Q)/\   positions  appropriate  for  CaMg(C03)2  suggests  that  there  are  basal  sequences 
of  limited  extent  within  the  crystals  that  have  essentially  the  ideal  structure  and 
composition  appropriate  to  the  dolomite  formula.    The  various  models  that  we  are 
about  to  test  afford  a  means  of  setting  aside  a  portion  of  the  total  carbonate  pack- 
ets as  an  ideal  alternation  CO3— Ca— CO3— Mg-C03  (henceforth  referred  to  as  the 
"ideally  alternating  portion"),  of  arranging  the  remaining  packets  (henceforth  refer- 
red to  as  the  "residuum")  in  various  ways,  and  then  of  mixing  the  two  portions  to 
various  degrees.    We  restrict  ourselves  to  one  bulk  composition,   55  mol  percent 
CaC03. 

Possible  ways  of  distributing  excess  CaC03  in  the  residium  include:  (1) 
single  or  multiple  excess  Ca  planes  inserted  at  regular  intervals  in  an  ideally  al- 
ternating cation-plane  succession  (regular  mixed  layering);  (2)  single  or  multiple 
excess  Ca  planes  inserted  at  random  in  an  ideally  alternating  cation-plane  succes- 
sion (one  kind  of  random  mixed  layering);  (3)  single  or  multiple  excess  Ca  planes 
inserted  at  random  in  a  random  succession  of  cation  planes  (another  kind  of  random 
mixed  layering);  (4)  solid  solution  in  the  Mg  cation  planes  of  an  ideally  alternating 
cation-plane  succession;  (5)  solid  solution  in  the  Mg  cation  planes  of  a  random  sue- 
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cession  of  cation  planes.     In  subsequent  discussion  we  will  refer,   for  example,   to 
"Type  3"  residua  and  "Category  3"  (of  mixed  layering). 

The  last  three  types  have  already  been  treated  by  the  two-packet  model,  ex- 
cept that  there  was  no  ideally  alternating  portion  to  be  considered  in  addition,  and 
that  for  Type  3  residua,  one  could  not  insert  multiple  excess  Ca  planes  at  random, 
i.e.  set  a  lower  limit  to  the  length  of  an  inserted  Ca-plane  sequence. 

The  repeated  reference  in  the  second  preceding  paragraph  to  "single  or  mul- 
tiple" reflects  our  interest  in  testing  the  effect  upon  the  diffraction  diagram  of  suc- 
cessively longer  series  of  neighboring  Ca-populated  cation  planes.     The  extreme 
reluctance  even  at  832°  C  of  Ca-rich  sedimentary  dolomite  crystals  to  reorganize 
to  an  equilibrium  two-phase  assemblage  (Goldsmith  and  Graf,    1958)  suggests  that 
the  ideal  basal-plane  succession  is  closely  approximated  and  that  long  series  of 
Ca  planes  are  highly  improbable.    In  general,  we  attempt  to  prevent  Mg-populated 
planes  from  being  nearest  cation-plane  neighbors  because  most  of  the  natural  do- 
lomites examined  thus  far  that  deviate  from  ideal  composition  are  Ca-rich,  and  be- 
cause we  anticipate  that  there  will  not  be  enough  variable  parameters  in  the  three - 
packet  models  to  allow  for  both  multiple  Ca-plane  sequences  and  multiple  Mg-plane 
sequences. 

Some  Three-Packet  Markov- 1  Models 

Table  1  gives  packet  types  for  four  three-packet  models   that  differ  only 
in  the  thicknesses  of  the  Type  3  packets  used  to  introduce  into  the  crystal  the  excess 

TABLE   1   -  PACKET  TYPE  DEFINITIONS   FOR  SOME  THREE-PACKET  MODELS^ 


Model 

Packet  type 

A 

B 

C 

D 

CO3 

CO3 

CO3 

CO3 

1 

Ca 

Ca 

Ca 

Ca 

CO3 

CO3 

CO3 

CO3 

Mg 

Mg 

Mg 

Mg 

CO3 

CO3 

CO3 

CO3 

2 

Mg 

Mg 

Mg 

Mg 

CO3 

CO3 

CO3 

CO3 

Ca 

Ca 

Ca 

Ca 

CO3 

CO3 

CO3 

CO3 

Ca 

Ca 

Ca 

Ca 

CO3 

CO3 

CO3 

3 

Ca 

Ca 

Ca 

CO3 

CO3 

Ca 

Ca 
CO3 

Ca 

These  models  are  devices   for  obtaining  Markov-2  dependence  using  Markov-1  pro- 
grams.    A  formal  development  of  this  relationship   is  given  in  Blyth,   Graf, 
and  Stemmler    (see  references) .     Model  B,    in  particular,    is  a   special  case  of 
the   four-packet  Markov-1    (two-packet  Markov-2)   model  of  Blyth,   Graf,   and 
Stemmler  with  packet  type  4,   CO3 — Mg — CO3 — Mg,    suppressed. 
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5  mol  percent  CaC03  above  50  mol  percent.     Packets  of  Types  1  and  2  are  identical 
in  composition  but  have  opposed  cation-plane  successions. 

Consider  Model  B  as  an  example.    The  expressions  for  amplitudes  scattered 
from  these  thicker  packets  have  correspondingly  more  terms,  e.g.  for  a  Type  2  pack- 
et,  summed  at  the  uppermost  carbonate  layer: 

V^^^  =  f^Q   (sin^)/X    |l  +  exp    [-STria^Csin  ^  )/x]  j    +    f^   (sin^)/X 

|exp    [-47ria2(sin^)/x]  j    +f^g(sin^)/X 
jexp    r-477-i(2a  +a  )(sin^  )/X  1  I  . 

There  are  now  nine  phase  shifts,  (p^^^l,  to  be  considered,    some  of  which 
are  identical  for  Model  B: 

^(11)^^(21)^^(31)^  ^(12)^^(22)^  ^(32)^g^(^^^^.^^j/^ 
^(13)^^(23)^^(33)^^g^^^(^.^^)/^^ 


The  Q^  matrix  becomes 

^(1)        f   ■  ,(ll)x  r(2)        ,  .  ,(12),          ,(3)        ,    .   ,(13), \        / 

c^^f'  'exp(-i(^      ')  0^2^      exp(-i,^'     ')  c^^f'  'exp(-i<^'     ')  \         q^^  q^^  q^^ 

o       I            f(l)        I    -J^Dv  ,(2)         ,    .   ,(22).           ,(3)        ,    .  ,(23). 

Qj=     °21        exp(-i0        )  C22f      exp(-i0        )  C23f      exp{-icp        ) 


c(l) 


,(31) 


:(2) 


,(32) 


(3), 


.  ,(33) 


c^f      exp(-ic^        )   C32f      exp(-i<^        )   Cggf      exp(-i<^        )/        \a^^  q^^  q/. 

where  the  Cg^  as  before  are  multipliers  showing  how  the  probability  f^^  is  changed 
by  the  knowledge  that  the  previous  packet  is  of  type  s.  The  eigenvalues  of  Qi  are 
obtained  by  calculating  the  coefficients  of  the  characteristic  polynomial 


=  1    ^21   ^22   ''23 


-X^  ^   ('5ll^^22^^33^^^  ^   (^21^12^S2^23^^31^13-^11^22-^llS3-^22^33^^ 

■'^ll'^22^33^^12^23'^31^''l3Sl''32-^12^21^33"^11^32^23"^22''l3^3r 

There  are  now  nine  equations  analogous  to  (11)  and  in  each  set  of  three  equations 
one  a^j  can  be  chosen  arbitrarily  (but  not  zero). 

The  three -packet  equivalent  of  equation  (13)  is  strictly  analogous: 


3,,,,, ,2       3K(R^^^)q^^V|q(^^|   W^^^ 
s=l  s=l  |l-Q      I 


In  the  three-packet  analogue  of  (5),  the  conditions  that  the  row  sums  be 
each  1  and  the  probabilities  in  the  crystal  be  stationary  become 
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^31^ 


^13^ 


(1) 
(1) 
(1) 
(1) 
(1) 
(1) 


-.2^'"-3/^'- 


Subtracting  (17)  from  (14)  and  rearranging  gives 

(c 
Similarly,   for  (18)  and  (15), 

(c 


./^(3)     .  ./Jl) 


and  for  (19)  and  (16), 


(^31-13)/^^'^=  (^23-32)/^^'^' 


(14) 
(15) 
(16) 
(17) 
(18) 
(19) 


giving 


(c 


12   "21 


'/f"'=fc23-32'/f'"-'=31-13'/f'"  =  ^- 


(20) 


If  the  packet  succession  is  random,  all  Cj^j's  have  the  value  one  and  k  would  be  ze- 

case.    All  the  c^.-'s  can  be 


ro,  although  the  derivation  is  not  valid  for  that  limiting 

,.  k.  f(i),  f(2\  3^d  f(3r 


given  in  terms  of  C|i,   C22»  c^g, 

Finally,  we  define  a  run  as  the  number  of  successive  cation  planes  of  a  giv- 
en kind  lying  between  two  cation  planes  of  the  other  kind  and  we  set  forth  in  table 
2  the  runs  produced  by  different  packet  successions  in  Models  A,  B,  C,  and  D.    Our 
manner  of  defining  Type    1  and  2  packets  has  the  useful  consequence  that  Mg  can- 
not occur  in  runs  longer  than  2. 


Preliminary  Calculations 

We  have  obtained,   by  hand  calculation  of  one -dimensional  Fourier  series, 
the  basal  spectra  for  a  number  of  Type  1  residua  (fig.  7).    These  intensities  as  shown 
do  not  constitute  basal  diffraction  diagrams.    Each  intensity  in  a  given  diagram  should 
be  reduced  (ignoring  mass  absorption  and  other  corrections)  to  a  percentage  "x"  of 
itself  given  by  the  figure  in  the  upper  right  hand  corner,  and  the  four  dolomite  reflec- 
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tions  should  be  drawn  in  at  intensities  (1-x)  percent  of  those  shown  for  dolomite. 
The  (00-3)  dolomite  reflection,  visible  on  most  but  not  all  of  the  single -crystal  di- 
agrams of  Goldsmith  and  Graf  (1958),   serves  as  a  convenient  intensity  cutoff.     Omit- 
ted from  figure  7  are  those  calculated  intensities  that  fall  within  log  units  below 
that  in  which  the  calculated  value  for  the  dolomite   |00-3J  reflection  lies. 

To  the  right  in  figure  7,  the  basal  diffraction  diagram  of  calcite  is  approached 
as  a  limit;  downward,  that  of  dolomite;  diagonally  downward  to  the  right,  that  of  a 
crystal  composed  of  a  calcite  slab  and  a  dolomite  slab  with  the  same  e-axis.    The 
approach  to  these  limiting  diagrams  as  a  function  of  repeating-sequence  length  is, 
however,  quite  slow.     The  computed  diagrams  of  figure  7  all  show  additional  reflec- 
tions, lying  between  the  four  ranges  defined  by  the  short  horizontal  bars,  that  are 
more  intense  than  the  dolomite    |00-3J   reflection,  even  after  making  a  correction 
for  the  "x"  factor.     If  the  array  of  diagrams  were  continued  far  enough  to  the  right 
or  diagonally  downward  to  the  right,  diffracted  intensity  between  the  four  defined 
basal  reflection  regions  would  eventually  sink  below  the  cutoff  value,  but  long  be- 
fore this  point  was  reached  there  would  be  strong  maxima  at  essentially  the   |00-6| 
and    JOO'IZJ    calcite  positions,   as  there  already  are  for  CasCCaMg)  and  Ca^lCaMgl^. 

The  simple  diagrams  of  figure  7,  therefore,  do  not  correspond  to  the  diffrac- 
tion features  we  seek  to  explain.     However,   if  we    were  to  set  up  longer  and  longer 
sequences  in  which  the  excess  Ca  planes  were  inserted  in  various  multiples  and/or 
at  irregular  intervals,  e.g.  Ca(CaMg)„Ca-(CaMg)Ca^(CaMg)r ,  we  would  approach 
a  Type  2  residuum  as  a  limit.    As  the  residuum  was  allowed  to  become  a  larger  and 
larger  percentage  of  the  total  crystal,  extra  reflections  other  than  those  of  dolomite 
would  progressively  attenuate  until  they  disappeared  into  the  background,  just  as 
the    JOO-SJ  and   |00-9J  reflections  did  in  the  two-packet  calculations.     Here,  how- 
ever,  because  the  (CaMg)  units  have  been  maintained  intact,   the  residuum  at  the 
limit  would  still  be  an  alternating  cation-plane  sequence,  although  an  imperfect 
one.    The  residuum,  therefore,  would  make  a  coherent  but  somewhat  weakened  con- 
tribution near  the  dolomite  position  in  all  four  ranges.     However,   calculations  for 
these  long  discrete  sequences  become  very  time  consuming,  and  we  are  forced  to 
Hendricks -Teller  three-packet  models  in  which  the  nature  of  the  sequence  is  de- 
fined by  statistical  parameters. 

The  value  of  calculations  like  those  upon  which  figure  7  is  based  is  further 
limited  by  the  absence  of  any  means  of  mixing  residua  and  ideally  alternating  por- 
tions.    Even  though  we  have  no  quantitative  knowledge  of  the  extent  of  that  mixing, 
we  know  that  it  must  occur  if  reflections  with  no  basal  component  are  not  to  be  bi- 
modal.    Such  mixing  will,   of  course,   be  an  additional  factor  weakening  the  extra 
residuum  reflections.    Information  of  the  type  contained  in  figure  7  did  prove  use- 
ful in  testing  computer  programs  and  in  recognizing  extra  residuum  reflections  that 
persist  in  Hendricks -Teller  curves  for  certain  parameter  values. 

Three-Packet  Markov- 1  Calculations 

Within  each  model  (A,   B,  C,  D  of  table  1)   the  relative  magnitudes  of  f(l) 
and  f'^'  can  be  varied,   subject  to  the  restriction  that  f'-^'+   f'^'  remains  constant, 
without  changing  the  bulk  composition.    We  choose  the  c^j's  so  that  Type  1  and 
Type  3  packets  will  mix  but  those  of  Type  2  will  be  largely  segregated.    Varying  the 
magnitude  of  f" '  thus  varies  the  percentage  of  the  crystal  into  which  a  fixed  Ca  ex- 
cess is  mixed  and  corresponds  to  moving  along  the  appropriate  column  of  figure  7. 
The  ci2  parameter  measures  the  degree  of  mixing  of  the  residuum  with  the  ideal  do- 
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oMg  (30)       COgCoMg  (20)      COjCoMq       (ISZ/j)      Co^CaMq  (15)      COjCoMq  (14)  Colcite,  Co  (10) 


Ca(CoMq)  (50) 


Ca(CoMq)j         (70) 


Ca(CaMq) 


COj(CoMq)2         (30) 


Dolomite,  CoMq 

1 

Figure  7  -  Basal  reflection  positions  and  intensities  computed  for  residua  having 
the  indicated  cation-plane  successions.     The  expression  Ca2(CaMg)3  repre- 
sents,  for  example,   the  regularly  repeating  succession  CaCaCaMgCaMgCaMg, 
The  intensities,   which  are  plotted  on  a  5-cycle  logarithmic  scale,   have  been 
normalized  to  a  constant  total  number  of  cation  planes  and  include  Lorentz  and 
polarization  corrections  appropriate  for  Cr  radiation  and  zero-layer  Weissen- 
berg  reflections.    The  four  short  horizontal  bars  repeated  at  the  base  of  each 
block  of  the  diagram  are  the  (sin5)/X  ranges  at  which  the    JOO-SJ  ,    J00-6|, 
|00-9|  ,   and  l00-12|   reflections  occur  for  the  compositional  range  calcite- 
magnesite.    The  number  In  parentheses  in  each  block  is  the  percent  of  cation 
planes  in  the  total  crystal  that  is  in  the  residuum  for  that  specified  residuum 
cation-plane  succession,   for  a  bulk  composition  of  55  mol  percent  CaC03  and 
for  cation  planes  without  solid  solution. 
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lomite  portion;  we  choose  to  make  calculations  at  Ci2=  .001  and  .1.    The  former 
value  is  a  minimal  one  that  should  correspond  to  the  two-slab  crystal  situation  of 
figure  7.    The  latter,   still  small  in  comparison  with  maximum  possible  0^2' s,   is  a 
first  guess  at  a  value  for  which  reflections  with  no  basal  component  might  become 
symmetrical,  yet  basal  reflections  still  be  asymmetrical  or  doubled. 

We  wish  our  c^j's  to  be  such  that  Ca  in   Model  A  occurs  only  in  runs  of  2 
(and  1),   in  Model  B  in  runs  of  3  (and  1),   in  Model  C  in  runs  of  4  (and  1),  and  in 
Model  D  in  runs  of  5  (and  1).     Furthermore,  we  would  like  to  place  Mg  in  runs  of  1. 

These  several  conditions  on  the  choice  of  c^j's  can  be  realized  in  all  four 
models  by  maximizing  C22>  ^^3,  and  C3P   by  minimizing  C21  (to  place  Mg  in  runs 
of  1)  and  C33  (to  forbid  Ca  runs  that  are  long  and  of  uncertain  length),  and  by  test- 
ing C]^2  ^t  '^01  ^^^  •^'    Packet  successions  131  thereby  generate  the  desired  runs 
in  Ca  (other  than  1).     Because  we  have  placed  restrictions  on  six  c^j's  (whereas 
only  four  can  be  chosen  independently  if  the  f  s  are  fixed),  and  because  the  allowed 
c^j's  are  functions  of  the  f  s,    a  certain  amount  of  preliminary  juggling  is  required 
to  get  the  best  compromise.     In  practice,  the  whole  set  of  Cj^j'  s  is  computed  to  be 
certain  that  none  are  negative  or  zero  (zero  c^^^'s  appearing  in  equation  (12)  or  the 
analogous  three-packet  expressions  would  lead  to  indeterminate  a^j's).    These  sets 
of  parameters  for  our  computed  three-packet  curves  are  given  in  table  3.    The  de- 
sired c^j  combinations  appear  to  have  been  satisfactorily  approximated,  except  for 
some  C2  j'  s  approaching  .  1,   an  order  of  magnitude  above  what  could  be  considered 
minimal.    The  effect  of  the  relative  m&gnitudes  of  the  f'  s  on  the  attainable  c^j*  s 
appears,  for  example,   in  Curve  32  where  one  is  disturbed  on  first  inspection  oy  the 
relatively  large  cj^j  and  the  relatively  small  0^3 . 

Setting  each  of  the  three  quotients  in  equation  (20)  in  turn  equal  to  k,   and 
substituting  maximum  and  minimum  values  (l/f'^',   0)  of  the  appropriate  Cjj's  for  a 
given  column  of  table  3,  three  apparent  ranges  are  obtained  (in  general,  different) 
over  which  k  may  extend.     The  real  allowed  range,  quoted  at  the  bottom  of  each 
column  in  the  table,   is  that  common  to  all  three  of  the  computed  ranges.     In  gen- 
eral, the  values  of  k  actually  used  in  computation  lie  close  to  zero  as  a  consequence 
of  making  C3J  and  0^3  large,  C23  and  C32  small,  and  cj2  and  C21  small.     Thus, 
^31~^13»  ^23~^32'  3ndci2-C2i,  but  this  relation  is  not  a  requirement  of  the  mod- 
el, which  can  treat  the  range  of  permissible  k's. 

An  increase  of  0^2  to  .1  for  curves  2,   10,    18,  and  26,  which  have  low  val- 
ues of  f^^',  was  achieved  by  allowing  k  to  increase  moderately.    The  somewhat  ov- 
erlarge  C2i's  referred  to  above  could  have  been  prevented  by  increasing  k,  which 
would  have  decreased  C21  but  maintained  the  necessary  transitions  into  residuum 
progressions  through  a  larger  C23.    Larger  k's  correspond,  of  course,  to  crystals 
progressively  more  polar  in  cation  distribution,  a  feature  that  might  be  capable  of 
experimental  verification  through  piezoelectric  tests. 

The  curves  are  reproduced  in  figures  8-11.    Also  included  are  curves  for  suc- 
cessions in  which  Type  3  packets  are  randomly  distributed  throughout  the  residuum. 
These  were  generated  by  modifying  the  parameters  of  Curves  2,  4,   6,  and  8  to  yield 
the  columns  of  parameters  denoted  by  2A,  4A,  6A,  and  8A  in  table  3.    Mixing  be- 
tween residuum  and  ideally  alternating  portion  has  been  held  at  the  same  level,  i.e. 
^12~  •^'  ^22  ^^  ^*^^^  maximized  at  the  expense  of  C21  and  C23.  and  C32  is  still  min- 
imized.   The  two  parameters  in  each  of  the  pairs  c^^,  Cj^3  and  03^,  C33  have  been 
brought  to  identical  values;  if  there  were  not  restrictions  remaining  upon  Cj2  ^i^d 
C32,  these  identical  values  would  be  1.    Because  C33  is  no  longer  trivial,  Ca  runs 
of  all  lengths  in  residua  are  possible  for  these  curves  based  upon  Model  A.    The 
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analogous  successions  based  upon  Model  B  would  have  Ca  run  lengths  of  1,   3,   5, 

7,    .  .  .  allowed;  Model  C,    1,   4,   7,    10,    ...    ;  Model  D,    1,   5,   9,    13,    We 

have  also  calculated  curves  26A,   28A,   30A,   and  32A,   based  upon  Model  D, 

The  curves  with  c;[2~  .001  reproduce  the  positions  of  the  maxima  of  figure 
7  with  high  precision.     This  can  be  seen  on  curves  such  as  1  and  9  where  the  f  s 
chosen  yield  residua  that  duplicate  those  of  CaCaMg  and  Ca2CaMg,  respectively. 
A  few  extra  residuum  reflections  occur  on  the  parts  of  these  curves  shown;  the  re- 
laxation of  On  2  to  .1   in  Curves  2  and  10  weakens  but  does  not  destroy  them.    As 
expected,   they  are  destroyed  in  curves  with  Type  2  residua  (compare  Curves  26  and 
26A).    At  very  large  f"'  values,  there  is  little  or  no  difference  between  the  curves 
for  crystals  with  Type  1  and  Type  2  residua  (compare  Curves  8  and  8A,   32  and  32A), 
because  the  Type  3  packets  are  so  outnumbered  by  Type  1  packets  that  it  makes 
very  little  difference  whether  they  are  dispersed  randomly  or  at  regular  intervals. 

The  extreme  sharpness  of  many  of  the  maxima  arising  from  ideally  alternat- 
ing portions  is  a  consequence  of  our  having  directed  parameter  choices  toward  de- 
fining the  residuum.    As  a  result,  the  principal  factor  preventing  these  sharp  maxi- 
ma from  having  infinite  intensity  is  the  nonzero  nature  of  C22>     Obviously,   it  would 
have  been  possible  to  define  models  giving  diffuse  maxima  at  the  dolomite  position 
and  sharp  residuum  peaks.     It  also  should  be  noted  that  models  such  as  A  and  C, 
in  which  packets  of  the  third  type  contain  an  odd  number  of  Ca  planes,  generate 
|00-9|  profiles  having  the  residuum  and  ideally  alternating  components  reversed 
in  position  from  that  in  I  00- el    and    |00-12|. 

The  largest  possible  component  separation  in  bimodal    |00-  12j    maxima,  that 
between  the  positions  for  calcite  and  dolomite,   is  realizable  (figs.  8-11)  for  an  in- 
finite range  of  Type  3  packet  thicknesses.     On  the  other  hand,  at  the  c-]^2  value  of 
.1,  which  we  have  been  using  as  reasonable  for  natural  crystals,  the  smallest 
possible  separation  decreases  as  the  Type  3  packet  thickness  decreases. 

We  have  not  made  calculations  for  our  three -packet  models  with  substitu- 
tion of  one  cation  into  basal  planes  occupied  by  the  other.    Given  the  amplitude 
contributions  to  structure  factors  for  dolomite  (Graf,    1961)  reproduced  here  in  table 
4,  together  with  their  numerical  equivalents  for  Cu  radiation,  we  can  predict  that 
for  a  fixed  i^^',   substitution  of  the  smaller,   more  weakly  scattering  Mg  atom  into 
Ca  planes  would: 

(1)  Shift  intensity  maxima  of  the  ideally  alternating  components  of  all  orders 
of  basal  reflection  to  higher  (sin  9)/\  values; 

(2)  Decrease  the  intensity  of  this  component  for   |00-9J    and   |00-12|,   but 
increase  it  for    |00*6|'  ; 

(3)  Decrease  the  intensity  of  the  ideally  alternating  component  of    |00-3|    up 
to  about  7  atomic  percent  Mg  substitution,   but  increase  it  thereafter. 

At  fixed  bulk  composition,   f^^',   and  Type  3  packet  definition,   substitution  of  Mg 
into  Ca  planes  would: 

(1)    Shift  intensity  maxima  of  residuum  components  of  all  orders  of  basal 
reflection  to  lower  (sin^)/x  values.    Although  both  Type  1  and  Type  3 
packets  are  more  Mg-rich  than  before,   fv3)  has  to  be  increased  to  main- 
tain the  bulk  composition.    Viewed  another  way,   the  ideally  alternating 
portion  withdraws  proportionately  more  Mg  relative  to  Ca  from  the  same 
bulk  composition  than  it  did  before  substitution,  leaving  a  relatively 
more  Ca-rich  residuum.     Thus,   Mg  substitution  into  Ca  planes,   in  the 
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TABLE  4  -  AMPLITUDE  CONTRIBUTIONS  TO  STRUCTURE  FACTORS 
FOR  THE  FIRST  FOUR  ORDERS  OF  DOLOMITE  BASAL  REFLECTION^ 


Basal 

reflection 

Amplitude   contribution 

■  00-3} 

f_     -   f        -  0.2444   f„   -  0.6771    f_ 
Ca           Mg                           C                           0 

17.56   -   10.84  -   1.28   -  4.97=  +0.47 

•  00- 6 j 

J 

f„     +    f        -   1.9402    f„   -   5.8472    f^ 
Ca          Mg                          C                          0 

14.67+     9.25   -   7.37   -  34.26  =  -17.71 

{00-9} 

f„      -   f        +   0.7187    f„   +   1.9969    f^ 
Ca          Mg                           C                           0 

12.18   -   7.78  +    1.95  +    8.69=  +15.04 

•  00-12  ■ 

f„     +   f„     +   1.7646    f„  +    5.3964   f 
Ca          Mg                           C                           0 

10.09+    6.53+    3.60+    17.43=  +37.65 

■'■  (After  Graf,  1961).   The  numerical  equivalents  are 

for  Cu  radiation  and  the  assumption  of  neutral  atoms. 


situation  we  have  defined,  would  lead  to  increased  separation  of  the 
two  components  of  each  reflection; 

(2)  Increase  the  intensity  of  the  residuum  component  for    |00  -12 j  ,   but  de- 
crease it  for  |00-6J  .     The  residuum  contains  the  same  number  of  cation 
planes  as  before,   and  the  Ca  planes  are  now  Mg- substituted,   but  they 
are  a  greater  fraction  of  the  total  cation  planes  in  the  residuum  so  that 
the  residuum  composition  is  more  Ca-rich  than  before; 

(3)  Decrease  the  intensity  of  the  residuum  component  for  |00-9|   and  in- 
crease it  for    |00-3|  ,     Substitution  of  Mg  weakens  scattering  from  Ca 
planes  in  Ca-plane  sites;  the  increased  number  of  Type  3  packets  pla- 
ces some  Ca-rich  planes  in  Mg-plane  sites  and  thus  increases  the  av- 
erage scattering  from  Mg-plane  sites. 

Analogous  reasoning  for  the  substitution  of  Ca  into  Mg  planes  shows  that 
intensity  changes  for  both  components  of    |00-3|    and    |00-9J   would  be  in  the   same 
direction  as  for  the  Mg-into-Ca-plane  case,   but  that  the  other  relations  discussed 
above  would  be  reversed.     In  particular,  the  two  components  of  each  basal  reflec- 
tion would  move  closer  together.     In  this  case,   for  fixed  f^^',   f'*^'  decreases  on 
substitution. 

Simultaneous  substitution  of  equal  atomic  percentages  of  Mg-into-Ca  planes 
and  Ca-into-Mg  planes  would,   for  our  assumptions,  decrease  the  intensity  of  both 
|00-9|   components,   increase  that  of  the    |00-3|   residuum  component,  decrease 
that  of  the    JOO  -31    ideally  alternating  component  up  to  about  3  atomic  percent  ca- 
tion substitution  and  then  increase  it,  and  make  no  significant  change  in  the  other 
quantities  discussed  above.     For  this  case,  with  fixed  f^  '  and  bulk  composition, 
the  ideally  alternating  portion  and  the  residuum  each  retain  their  composition  before 
substitution,   but  f^"''  must  decrease  and  f^  '  increase. 
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Some  Natural  Examples 

We  have  re-examined  the  group  of  zero-level  a-axis  Weissenberg  diagrams 
from  which  those  in  plate  2  of  Goldsmith  and  Graf  (1958)  were  chosen  for  reproduc- 
tion.    Because  these  dolomite  crystals  are  sharply  defined  replacive  rhombohedra 
from  sedimentary  limestones  of  Ordovician  and  Mississippian  age,  and  because  fer- 
ricyanide  staining  tests  show  that  they  retain  sharply  defined  compositional  zones, 
we  assume  that  any  basal  mixed  layering  observed  is  an  original  growth  feature  and 
does  not  involve  later  reorganization.     The  presence  of  diffracted  intensity  from 
Kci2  radiation,   out  to  the  iQ  region  of  (03-  0)  and  (00  -12),  does  not  noticeably  alter 
the  intensity  distribution  within  individual  components  of  the  doubled  reflections 
(although  it  broadens  each  component  somewhat),  and  we  have  ignored  this  factor 
in  making  visual  estimates.    We  note  in  general  that: 

(1)  There  are  no  basal  reflections  other  than  (00-3),   (00-6),   (00-9),    (00-12), 
and  (00-  15).    Although  the  presence  of  white  radiation  streaks  makes  it 
more  difficult  to  detect  weak  reflections,  we  conclude  that  any  of  the 
other  residuum  ordering  reflections   of  figure   7  that  are  at  least  as 
intense  as  (00-3)  must  be  considered  missing; 

(2)  The  lower  angle  components  of  basal  reflections,  when  distinct  com- 
ponents can  be  defined,   are  always  sharper  than  the  higher  angle  ones, 
and  sometimes  as  sharp  as  (03-0); 

(3)  Spacing  measurements,   calculated  for  Kct  radiation  without  film  shrink- 
age correction,   show  that  the  sin  Q/\  separation  of  the  two  components 
of  (00-12)  is  small  and  rather  constant  (.  0022— .  0033); 

(4)  The  ratio  of  component  intensities  within  (00  -9),  within  the  limits  of 
visual  observation,   is  the  same  as  that  for  (00-12); 

(5)  The  intensity  diffracted  in  higher  angle  components,  as  earlier  sup- 
posed, reaches  as  a  limit  the  positions  for  dolomite  of  ideal  compo- 
sition and  structure,   and  the  lower  angle  components  correspond  to 
dolomite  with  some  54  to  59  mol  percent  CaCOg. 

Mixed -layering  Categories  3  and  5  (p.  31),  which  randomize  the  cation- 
plane  succession  of  the  residuum,  are  ruled  out  by  the  presence  of  both  compon- 
ents of  (00-9)  on  the  Weissenberg  diagrams. 

Two  kinds  of  evidence  argue  against  regular  mixed  layering  (Category  1). 
The  first,  which  was  mentioned  earlier,   involves  the  computed  diagrams  of  figure 
7.     In  all  of  these  diagrams  there  are  extra  basal  reflections,   in  addition  to  those 
characteristic  for  dolomite,   much  stronger  than  the  (00-3)  reflection  detectable 
above  the  white  radiation  streaks  in  the  Weissenberg  diagram. 

The  second  argument  involves  a  consideration  of  which  reflection  compo- 
nents form  integral  series  of  orders.     For  convenience  in  comparing  Categories  1 
and  4,  we  assign  letters  (fig.   12)  to  the  basal  diffraction  features  involved  in  the 
Weissenberg  diagrams.     The  separation  of  (00-12)  components  in  these  diagrams  is 
so  small  that  we  can  immediately  focus  our  attention  on  the  curves  of  Model  A  (fig. 
8).     Of  those  curves  in  the  group  2,   4,   6,   8,    10,    12,    14,    16  that  show  a  doubled 
(OO  -121  ,  the  smallest  {sinQ)/\  separation  for  Model  A  is  .0034  (Curve  6)  and  for 
Model  B  is  .0054  (Curve  12).     Model  A  predicts  basal  reflections  in  which  the  ser- 
ies of  ideally  alternating  components  CDG  forms  a  rational  series  of  orders  (con- 
stant (sin  B)/\  separation),   and  the  residuum  components  BEF  do  not.     Thus,  the 
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D      E 


(00-3) 


(00'6) 


(00-9) 


(00-12) 


Figure  12  -  A  schematic  representation  of  the  basal  diffraction  features  observed 
In  a  group  of  zero-level  a-axis  Weissenberg  diagrams.  The  letters  represent 
the  {sin  9  )/\  positions  at  which  the  several  intensity  maxima  occur. 

separations  C-D  and  D-G  of  Curve  4,    Model  A,   are  equal  (table  5),   except  for  a 
small  effect  from  c-io  being  .  1  rather  than  zero.     On  the  other  hand,  a  residuum  with 
excess  Ca  in  solid  solution  in  Mg  basal  planes  and  with  ideal  cation-plane  alter- 
nation would  have  BDF  as  a  rational  series    (assuming  only  moderate  mixing  between 
residuum  and  ideally  alternating  portion)  and  the  higher  angle  components  CEG  would 
form  another  rational  series  with  a  slightly  different  separation.     If  the  natural  crys- 
tals are  described  by  Category  1,  therefore,  the  entries  for  a  given  sample  in  the 
last  three  columns  of  table  5  ought  to  be  essentially  identical.     If  Category  4  holds, 
the  entries  in  the  first  three  columns  ought  to  be  essentially  the  same,   and  those  in 
the  middle  three  columns  ought  to  be  essentially  identical  but  slightly  larger  than 
those  in  the  first  three  columns.    A  complete  set  of  measurements  for  one  sample 
and  scattered  values  for  three  others  are  all  clearly  consistent  with  the  second  hy- 
pothesis, which  involves  Category  4. 

The  choice  between  Categories  2  and  4  remains  to  be  considered.     In  theory. 
Category  2  should  be  tested  by  model  calculations  in  which  c;^  ^  =  033  =  C;^3  =  Coi  =  1. 
However,  as  we  have  seen,  restrictions  upon  the  other  c^^j's  make  it  possible  in  the 
calculation  of  Curves  2A,   4A,   6A,   and  8A  to  reduce  these  four  parameters  only  to 
the  values  given  in  table  3.     Even  at  these  values,  and  with  0^2  set  at  the  low  fig- 
ure of  .1  compared  with  maximum  possible  c-|2's  of  1 . 5  to  5.5,  the    |00-9J   intensi- 
ty contributions  from  residuum  and  ideally  alternating  portion  in  all  four  curves  have 
merged  into  a  feature  having  only  a  single  maximum.     Mixed  layering  of  Category  2 
appears  inadequate  to  explain  single  crystals  with  bimodal  (00-9)  reflections,  al- 
though it  is  one  possible  explanation  for  those  with  asymmetrical  (00-9)  profiles. 

The  structural  model  for  replacive  sedimentary  dolomite  rhombohedra  that 
emerges  from  this  lengthy  argument  is  one  in  which  excess  Ca  is  accommodated  by 
substitution  in  Mg  layers,   and  regions  of  the  crystal  containing  essentially  pure 
Mg  layers  and  substituted  Mg  layers  alternate  on  a  fine  enough  scale  so  that  bas- 
al diffraction  features  characteristic  of  mixed  layering  are  observed.     This  conclu- 
sion is  in  agreement  with  our  earlier  observation  that  for  a  Ca-rich  bulk  composition 
Ca  substitution  into  Mg  planes  decreases  the  separation  of  components  of  bimodal 
basal  reflections.    It  follows  that  these  crystals  must  have  grown  under  conditions 
near  the  limit  of  the  field  within  which  dolomite  of  ideal  structure  and  composition 
can  form,   so  that  minor  fluctuations  of  environmental  variables  were  adequate  to 
move  the  system  back  and  forth  across  that  boundary.     The  incorporation  of  Ca  into 
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Mg  layers  of  dolomite  is  presumably  an  additional  manifestation  of  the  difficulty 
of  precipitating  anhydrous  MgCOg  in  sedimentary  environments. 

Our  small  group  of  three-packet  curves  does  not  explicitly  consider  impure 
cation  layers,   multiple  Mg-layer  runs,   and  poorly  alternating  sequences  of  Ca5oMg5o 
molar  composition,   and  the  physical  significance  of  fixing  c^2  at  .  1  is  not  well  un- 
derstood.    It  is  possible,   therefore,  that  there  are  other  favorable  structural  arrange- 
ments that  we  have  failed  to  consider,   and  the  model  just  described  must  be  regard- 
ed as  only  suggestive.    At  best,   it  applies  only  to  samples  giving  diffraction  features 
like  those  of  plate  2-C  in  Goldsmith  and  Graf  (1958),   a  sharp  (30-0)  reflection  and 
bimodal  (00-12)  and  (00-9)  maxima.     If  (00-9)  were  merely  asymmetrical  rather  than 
bimodal,   it  would  not  be  possible  with  the  kind  of  reasoning  we  have  used  here  to 
make  even  qualitative  statements  about  the  point  versus  layer  distribution  of  excess 
Ca.    If,   in  addition,    (30-0)  were  smeared  out,  the  possibility  of  the  rhombohedron 
being  a  mixture  of  discrete  phases  (macrozoning)  would  come  into  question. 


DISCUSSION 

The  geometric  properties  of  the  various  kinds  of  diagrams  that  can  be  used 
to  present  the  results  of  Hendricks -Teller  calculations  should  be  of  interest  in  the 
study  of  any  petrologic  system  in  which  rock-forming  minerals  show  interlayering 
and  solid  solution.     The  compounds  involved  in  the  interlayering  need  not  belong  to 
the  same  space  group,   but  should  be  structurally  related  to  the  extent  that  some  re- 
flections persist  through  the  whole  compositional  range.     Thus,  we  have  considered 
mixtures  of  basal  slices  of  the  calcite  and  dolomite  structures  and  could  have  in- 
cluded some  of  huntite  if  there  had  been  any  good  reason.    We  have  not  presented 


TABLE  5  -  VALUES  OF   A  (sin^)/X  CALCULATED  FROM  MEASUREMENTS  OF  FOUR 

ZERO-LEVEL  a-AXIS  WEISSENBERG   DIAGRAMS,    COMPARED  WITH  THE 

SAME  PARAMETERS   FROM  TWO  THREE-PACKET  MODELS^ 


(sin 

6)/\  separations  between  intensity  maxima 

Sample 
number 

A-B 

B-D 

D-F 

A-C     C-E     E-G 

A-C     C-D 

D-G 

D-19-8 

.0933 

.0933 

.0933 

.0943 

.0958 

D-23-10 

.0937 

.0936 

.0934 

.0938 

.0957 

D-21-12 

.0936 

.0934 

.0934 

.0955 

D-24-7 

.0928 

.0933 

.0931 

.0942    .0937    .0945 

.0942    .0919 

.0963 

Curve  4 
(fig.  8) 

.0957 

.0883 

.1123    .0746 

.0937 

.0934 

Curve  6 
(fig.  8) 

.0950 

.0898 

.0805 

.0928 

•*■  A  blank  results    from  one   of   the   lettered  quantities    either  being  unmeasurable 
or  not   readily   calculable  with   the   available   computer  programs.      The   signif- 
icance  of    the    letters    is   given   in   figure    12.      Sample   D-19-8    is    from   the  Mis- 
sissippian   Sedalia    ("Fern  Glen")    Dolomite,    Calhoun  County,    Illinois;    D-23-10 
is    from   the   Ordovician  Fort  Atkinson    (Divine)    Dolomite,    Maquoketa  Group, 
Divine,    Illinois;    D-21-12    is    from   the  Ordovician  Galena-Platteville  Dolomite, 
Oglesby,    Illinois;    and   D-24-7    is   a    cleavage   fragment   of   a   dolomite   from 
Guanajuato,    Mexico.    Because  no   film   shrinkage   correction  was   applied   to    the 
measurements,    their  absolute  values   cannot   be  used   to   estimate   Ca   content   of 
the  dolomites. 
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our  results  on  a  wide  variety  of  diagrams  because  there  are  not  yet  enough  experi- 
mentally measured  diffraction  profiles  of  the  pertinent  carbonates  for  us  to  be  sure 
that  such  diagrams  would  be  of  practical  use,  and  because  we  are  still  undecided 
about  how  to  treat  centroids.    There  are,  however,   several  points  about  the  prop- 
erties of  such  diagrams  that  should  be  mentioned: 

(1)  The  number  of  variable  parameters  that  can  usefully  be  defined  in  a  model 
is  limited  by  the  number  of  experimentally  measurable  quantities  against 
which  they  can  be  graphed.     The  optimum  experimental  method  would  be 
diffractometry  of  single  crystals  or  single-crystal  plates  to  yield  a  cen- 
troid  value  and  the  positions  of  either  one  or  two  maxima  for  every  reflec- 
tion,  including  ordering  reflections.    Parameters  involving  the  difference 
between  two  positions  in  the  same  reflection     (e.g.  centroid  and  intensity 
maximum,  or  two  intensity  maxima)  sacrifice  one  measurable  but  are  exper- 
imentally obtainable  without  elaborate  measuring  procedures  and  are  less 
sensitive  to  errors  in  assumed  packet  thickness  than  are  absolute  angular 
positions. 

(2)  The  several    |00-6|  diagrams  of  figure  4,   and  likewise  the  several    |00-12| 
diagrams,  are  quite  similar  geometrically,    A  rapid  but  approximate  method 
of  constructing  additional  diagrams  would  involve  suitable  end -point  shift- 
ing and  linear  compression  or  expansion  of  the  horizontal  and  vertical  scales 
on  one  of  the  existing  diagrams.    A  considerably  better  approximation  could 
be  made  if  the  new  diagram  lay  between  two  calculated  ones,   so  that  a  mod- 
ified geometry  could  be  sketched  by  interpolation.     These  methods  will  be 
most  successful  if  the  scattering  factors  of  the  several  pure  packets,   and 
hence  of  the  partially  cation-substituted  ones,  are  very  similar. 

(3)  It  might  be  possible  to  construct  diagrams  like  those  of  figure  4  more  rap- 
idly by  deriving  analytical  expressions  for  certain  geometrically  constrained 
points  within  the  diagrams  and  solving  for  these  points  in  a  particular  di- 
agram.   The  common  intersections  of  the  families  of  centroid  curves  and  max- 
imum-intensity-position curves  are  surely  not  geometrically  rigorous,   if 
only  because  the  several  atomic  scattering  factor  curves  are  not  precisely 
parallel  in  any  finite  range  of  (sin^  )/X  .     However,   for  systems  in  which 
packet  scattering  factors  are  no  more  dissimilar  than  in  the  CaMg  carbon- 
ates,  and  for  which  the  assumption  of  constant  packet  thicknesses  inde- 
pendent of  bulk  composition  is  acceptable,  common  intersection  is  an  ade- 
quate assumption.     In  practice,   this  would  mean  holding  scattering  factors 
constant,   independent  of  (sin^  )/X  ,   and  refining  them  and  the  solution  by 
successive  approximation  until  consistent. 

At  the  maximum- intensity-position  intersection  for  a  given  order  of  bas- 
al reflection,    (sin^  )/X  being  restricted  to  a  range  of  values  appropriate  for 
that  order,    (  d  lj\y/  d   (sin  9)/k)=  0  is  single-valued  in  f(l)  (or  f(2))if  c  is 
taken  to  be  either  0^3^  ^"^  ^"V  value  of  c  greater  than  that  c  below  which 
the  reflection  becomes  bimodal.     This  intersection  is  one  point  on  a  line 
that  is  the  locus  of  end  points  of  the  curves  in  the  upper  and  lower  fields  of 
a  given  diagram.    Throughout  the  rest  of  the  line,  the  expression  is  triple- 
valued  in  f(l)  (or  f(2))  for  a  given  c  below  that  c  at  which  the  reflection 
becomes  bimodal  (i,e,  there  is  a  minimum  lying  between  two  maxima,  and 
it  is  the  f '  s  of  the  maxima  that  are  sought) ,     The  smallest  value  of  c  at 
which  the  expression  is  still  single-valued  in  f(l'  (or  f(2))  is  that  of  the 
bounding  curve  of  the  envelope  of  continuous  maximum-intensity-position 
curves , 
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The  centroid  calculation  is  defined  by 
sin    ^n  ,  \        /  \       /Sin    ^- 


sin 


sin  9        sin  9\       sin    9\    I  X      /sin  9      sin  9^\      /sin 

9^  L  9^ 


X  X 

where  (sin  9  2)'^^   is  the  centroid  position  and  (sin  ^  ]^)/X  and  (sin  ^3)/Xare 
the  limits  between  which  intensity  is  measured.     Only  at  the  centroid  curve 
intersection  does  a  single  (sin^  2^^'^'   ^        (°'"  ^      )  satisfy  the  equation  for 
all  values  of  c.     The  chance  of  being  able  to  manipulate  a  complicated  in- 
tegral equation  is  much  poorer  than  for  the  partial  differential  equation  of 
the  preceding  paragraph.     Furthermore,    it  is  not  possible  to  secure  a  sim- 
plified solution  of  equation  (13)  for  c=  0,   as  would  be  desirable  for  defin- 
ing the  bounding  curve  of  the  envelope  of  centroid  curves,   because  the  Hen- 
dricks-Teller equation  is  not  valid  for  this  limiting  case.    A  completely  new 
expression  would  have  to  be  derived, 

(4)    To  use  diagrams  like  those  of  figure  4  for  estimating  solid  solution  in  two- 
packet  calculations,   one  has  to  interpolate  between  diagrams  for  both  of  these 
cation  substitution  variables.     If  one  bulk  composition  is  of  particular  im- 
portance,  sets  of  diagrams  at  that  fixed  composition  can  be  constructed, 
each  set  with  a  different  fixed  value  of  substitution  into  one  packet  type. 
The  parameter  c  and  the  degree  of  substitution  into  the  second  packet  type 
are  shown  as  ordinate  and  abscissa  on  individual  graphs,  which  can  then 
be  contoured  with  families  of  iso— (sin  9  )/X  curves  for  centroid  position, 
intensity  maximum  position,   separation  between  maxima  of  bimodal  peaks, 
centroid- maximum  separation,  etc.     With  these  graphs,   only  one  cation- 
substitution  variable  would  have  to  be  interpolated  between  diagrams. 

Evidence  of  carbonate  mixed  layering  obtained  from  measurements  not  in- 
volving X-ray  diffraction  would  be  particularly  valuable.     The  possibility  of  polar 
basal-plane  successions  could  be  tested  by  seeking  optical  activity,   pyroelectrici- 
ty,  and  appropriate  piezoelectric  moduli,   if  suitable  samples  could  be  prepared, 
Zen  (1962)  has  outlined  general  monitoring- solution  and  titration  techniques  for  test- 
ing whether  or  not  sufficiently  long  series  of  packets  of  the  same  kind  might  result 
in  polyphase  behavior. 
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APPENDIX 

Digital  Computer  Procedure 

ILLIAC  I,  which  was  taken  out  of  service  at  the  end  of  1962,  was  a  fixed- 
point  one-address  computer  having  a  high-speed  store  of  magnetic  cores  with  a 
capacity  of  1024  storage  locations   and  a  magnetic  drum  containing  12,800  addi- 
tional locations.     Each  storage  location  contained  40  binary  bits.     Information  was 
introduced  and  output  obtained  using  perforated  paper  tape. 

Each  computation  for  the  two-packet  model  took  approximately  1  second, 
and  each  for  the  three-packet  model  about  14  seconds.     The  accuracy  lost  by  fail- 
ing to  maintain  the  maximum  possible  number  of  significant  figures  (39  binary  bits, 
or  1 1  to  12  significant  figures)  at  all  stages  in  the  calculation  was  considered  ir- 
relevant because  of  the  small  number  of  significant  figures  to  which  atomic  scat- 
tering factors  are  known. 

Two-packet  calculations  utilized  both  preliminary  and  variable  internal 
scaling.     Three-packet  calculations  were  written  in  the  slower  floating  point  form 
because  it  was  necessary  at  one  place  to  obtain  the  eigenvalues  of  a  third  order 
matrix  with  complex  elements.    A  1.8th  order  iterative  procedure  was  used  to  find 
the  roots  of  the  characteristic  polynominal. 

In  a  typical  set  of  two-packet  calculations  in  which  it  was  desired  to  ob- 
tain   |00-6|   profiles  for  several  (f"',    f'^^,   c)  combinations,   a  series  of  100  (sin 
d)/\  values  at  an  interval  of  0,0005  was  read  into  the  drum,   together  with  associ- 
ated fcai   fjvlg'  ^'^'^  ^00-5  ■^slues.     (It  would  be  better  practice  in  a  revised  program 
to  store  a  table  of  this  information  at  some  selected  (sin  9)/\   interval,   and  require 
the  computer  to  carry  out  linear  interpolation.)    The  program  together  with  the  first 
(f'-"-',    f'2),  c)  was  read  into  the  high-speed  store  by  the  electronic  tape  reader  from 
a  master  tape.     The  first    (  (sin  0)/X,   fQg,   fy[g,   fcOo)  ^^^  then  transferred  to  the 
high-speed  store,  a  value  of  1  computed,   the  second  set  transferred,  etc.  ,   until 
all  100  points  had  been  calculated.     Then  a  new  (f'     ,   f'^',    c)  was  read  in  from  a 
data  tape,   and  the  whole  operation  repeated.     In  contrast,  three-packet  model  da- 
ta were  not  stored  on  the  drum,    so  that  after  a  change  of  parameters  (f'   ^    f'   ', 
f^     ,   cji,   C22>   ^33'   1^)  each     (  (sin  ^)/X,    f^^,   fj^ig'   ^03)    ^^^  obtained  in  turn 
from  tape. 

The  calculations  should  be  reprogrammed  for  FORTRAN  or  some  other  com- 
monly used  programming  system.  We  see  little  advantage  in  making  the  ILLIAC  I 
program  generally  available. 
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